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COMPUTER PROGRAM FOR THERMODYNAMIC ANALYSIS OF OPEN-CYCLE 
MULTISHAFT POWER SYSTEM WITH MULTIPLE REHEAT AND INTERCOOL 

by Arthur J. Classman 
Lewis Research Center 

SUMMARY 

This report presents a computer program for the thermodynamic analysis of an 
open-cycle multishaft power system with multiple intercooling and reheating capability. 
The program can compute specific power output, specific fuel consumption, and cycle 
efficiency for power systems having any number of shafts up to a maximum of five. On 
each shaft there can be any number of compressors and turbines up to a maximum of five 
each, along with any specified number of intervening intercoolers and reheaters. A re- 
cuperator can be included in the system. Turbine coolant flow can be accounted for. 

The combustion-gas thermodynamic properties are valid for any fuel consisting of hy- 
drogen and/or carbon only. The program should be used with maximum temperatures 
no higher than about 2000 K (3600° R) because molecular dissociation is not included in 
the stoichiometry. 

This report presents a description of the general open-cycle system that can be 
analyzed, the analysis method, and a complete description of the computer program. 

The description of the computer prc^ram includes a detailed explanation of input and 
output including an illustrative example, a discussion of the main program and each of 
the subprograms, the dictionary of program variables, and the program listii^. 


INTRODUCTION 

Open-cycle systems for groimd power-production applications, both stationary and 
propulsive, are currently being studied at the Lewis Research Center. One of the first 
phases in the analysis of any power cycle is the determination of the thermodynamic 
performance as a function of the cycle variables. Such cycle analyses are performed 
most readily by computer. Therefore, a computer program was written to perform the 
thermodynamic analysis of open-cycle power systems. 



The computer program can be used to analyze power systems having any number of 
shafts up to a maximum of five. On each shaft there can be any number of compressors 
arid turbines up to a maximum of five each, aloi^ with intervening intercoolers and re- 
heat combustors. A recuperator can be included in the system, and turbine coolant flow 
can be accounted for. Any fuel consisting entirely of hydrogen and/or carbon can be 
used. The program should be used with maximum temperatures no higher than about 
2000 K (3600° R) because molecular dissociation is not included in the stoichiometry. 

This report presents a description of the general open-cycle system that can be 
analyzed, the analysis method, and a complete description of the computer program. 

The computer program description includes a detailed explanation of input and output in- 
cluding a sample case, a discussion of the main program and each of the subprograms, 
a program listing, and a dictionary of program variables. 


SYSTEM ANALYSIS AND DESCRIPTION 

An example two-shaft power system is used to show the type of analysis that can be 
performed by the computer program being described in this report. The general open- 
cycle power system that can be analyzed is then described. 


Example-System Analysis 

An example of a power system that can be analyzed using the computer program de- 
scribed herein is the two-shaft open-cycle power system with intercooling and reheating 
that is shown schematically in figure 1 and thermodynamically on a temperature-entropy 
diagram in figure 2. Air enters the low-pressure compressor and is partially com- 
pressed. The air is then cooled within the intercooler. Since heat is generated during 
compression and compressor work is proportional to temperature, the use of intercool- 
ing serves to reduce the total compression work required for a given pressure ratio. 

The cooled air is then further compressed in the high-pressure compressor and heated 
in the combustor by the burning of fuel. Partial expansion occurs in the high-pressure 
turbine, which drives the high-pressure compressor, and heating again occurs, this 
time in the, reheater. Since turbine work is proportional to temperature, the use of re- 
heating increases the total work obtained for a given expansion. Final expansion to at- 
mospheric pressure then occurs in the low-pressure turbine, which drives both the low- 
pressure compressor and the generator. 

The computer program described in this report can be used to analyze such a power 
system to determine cycle efficiency as a function of turbine-inlet temperature and com- 
pressor pressure ratio. The improvements in cycle efficiency resulting from the use of 
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intercooling and reheating can also be determined. Cycle efficiency computed by the 
program is plotted against turbine-inlet temperature in figure 3 for a system with one 
intercooling and one reheating (fig. 1), a system with one intercooling only, and a simple 
system (no intercooling and no reheating). For these three systems an increase in 
turbine-inlet temperature from 1090 to 1370 K (1500° to 2000° F) increases cycle effi- 
ciency by about 20 to 25 percent; a further increase from 1370 to 1650 K (2000° to 
2500° F) increases cycle efficiency by another 10 to 15 percent. Operation at higher 
turbine-inlet temperatures, therefore, results in significant improvements in cycle effi- 
ciency. 

The improvements in cycle efficiency resulting from intercooling and reheating are 
also indicated in figure 3. At a temperature of 1370 K (2000° F) as an example, the use 
of one intercooUng alone increases the efficiency from 39 to 43 percent and the addition 
of one reheating along with the intercooling increases efficiency to 45 percent. For a 
plant of fixed power output, this corresponds to about a 10-percent reduction in fuel con- 
sumption with intercooling and about 5 -percent additional reduction in fuel consumption 
with the reheating added. At other temperatures, the benefits are about the same mag- 
nitude with efficiency reaching about 50 percent at 1650 K (2500° F). The reductions in 
fuel cost associated with such improvements in efficiency must be considered against the 
additional capital cost of the plant. Additional intercoolings and reheatings can be used 
to 3d eld further improvements in cycle efficiency. The incremental improvements, 
however, become smaller and smaller as additional intercool and reheat steps are 
added. The optimum number of intercoolers and reheaters to be used will depend on a 
detailed economic analysis for a given plant. 

In order to realize the benefits associated with the use of intercooling and reheating, 
the system must operate at greater pressure ratios than the simple cycle. This is indi- 
cated in figure 4, where cycle efficiencies for the simple system and for a system with 
one intercooling and one reheating are plotted against pressure ratio at a turbine-inlet 
temperature of 1370 K (2000° F), as an example. For the simple cycle the optimum 
pressure ratio is 30, and efficiency decreases markedly if the pressure ratio is less 
than 20 or greater than 40. For the cycle with intercooling and reheating, the optimum 
pressure ratio is about 120. However, the curve is very flat in the region of maximum 
efficiency, and it is possible to operate at a pressure ratio of. about 75 without a signifi- 
cant decrease in efficiency. This pressure ratio of 75 is still considerably higher than 
that for the simple cycle. Thus, with intercooling and reheating, not only will additional 
heat exchangers be required, but also additional stages will be required for the turbo- 
machinery. 
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General-System Description 


The example system just discussed is only one of many possible variations for open- 
cycle power systems. The computer program presented in this report can be used to 
analyze open-cycle power systems having different numbers and types of components. 
Figure 5(a) shows schematically a system having the maximum number of components 
that can be included in the analysis; figure 5(b) shows a system having the minimum 
number of components that must be included. Any system intermediate between these 
two can be analyzed. The size of the maximum system was selected on the basis of it 
being large enough to include all cases of possible interest. 

In figure 5(a) the system capability includes an inlet, a series of compressors and 
intercoolers, a recuperator, a primary burner, a series of turbines and reheaters, and 
a diffuser. There can be as many as five shafts. On each shaft there can be as many as 
five compressors with any specified number of intervening intercoolers and five turbines 
with any specified number of intervening reheaters. Output power can be obtained from 
each shaft. Part of the compressor exit flow can be used as turbine coolant or lost 
through recuperator leakage. The turbine coolant, which can be directed to any or all 
of the turbines, rejoins the primary flow at the exit of the turbine beii^ cooled. 

Any liquid or gaseous fuel consisting entirely of hydrogen and/or carbon can be 
used. A fuel boost pump or compressor, with intercooling if desired, can be provided 
with the required drive power taken from the system output. 

The subscripts shown for the compressors, intercoolers, burners, and turbines 
correspond to the double subscript notation used in the program. The first subscript 
refers to the component number, and the second subscript refers to the shaft number. 

An understanding of the proper order of subscripts is necessary to correctly prepare 
program input. Specific directions for properly specifying input-variable subscripts are 
given in the section Description of Input and Output. 


METHOD OF ANALYSIS 

The equations used for the open-cycle power system thermodynamic analysis are 
presented in this section. Equations are presented for the thermodynamic properties of 
the gas, analysis of the fluid state changes in each component, and computation of the 
system performance in terms of specific power output, specific fuel consumption, and 
cycle efficiency. The calculation l(^ic for combining the components into the desired 
system is discussed as needed. All symbols used in the analysis are defined in appen- 
dix A. 
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Thermodynamic Properties 


To perform a cycle thermodynamic analysis, it is necessary to relate enthalpy 
change and constant-pressure entropy change to the gas temperature change and com- 
position. This is done herein as follows: The heat capacity c of each component 
of the general combustion gas is expressed as a polynomial function of temperature. 
Then, from the stoichiometry associated with the reaction of 1 mass unit of air plus as- 
sociated inlet humidity with f mass units of fuel of composition CH^^, a single equation 
for the combustion-gas heat capacity is obtained. Finally, the enthalpy and constant- 
pressure entropy changes are obtained as 

/■T 

Ah= / ^c„dT (1) 



( 2 ) 


The equations for Ah and A(p as functions of initial T^ and final T2 tempera- 
tures, fuel-air ratio f, fuel composition (function of x), and absolute humidity of the 
inlet air are derived in appendix B. These equations are valid over a temperature range 
of 200 to 2000 K (360° to 3600° R). A maximum temperature of 2000 K (3600° R) was 
chosen because the combustion stoichiometry does not include molecular dissociation, 
which starts to appear at about this temperature. The basic assumption inherent in 
these thermodynamic equations, as well as in most of the analysis, is that the gas obeys 
the ideal gas law. 

For brevity in the analysis to follow, equation (B7) for Ah and equation (B12) for 
A<p will be denoted as 


Ah = H(Tj,T2,f) (3) 

and 

A<p = $(Tj,T2,f) (4) 

Since the fuel composition and inlet-air absolute hmnidity have sii^lar values for any 
given system analysis, they were not included as argiunents of these functions. 
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Components 


The equations used to analyze the flow through the inlet, compressors, intercoolers, 
burners, turbines, diffuser, recuperator, and fuel system are presented in this section. 

Inlet . - Since this is a ground-based power system, inlet total and static conditions 
were assumed equal. The inlet component would probably be some t3q)e of air filter. 
Therefore, for the flow leavir^ the inlet 


Pi = Po^i 

Wj = 1 + m 


(5) 

( 6 ) 
(7) 


The ambient static temperature Tq, ambient static pressure Pq, inlet pressure recov- 
ery rj, and inlet absolute humidity m are program inputs. 

Compressor . - The compressors and intercoolers are analyzed in flow sequence. 
The inlet conditions to the first compressor are those from the inlet component as speci- 
fied by equations (5) to (7). For any subsequent compressor, the inlet conditions are 
the exit conditions from the component, either a compressor or an intercooler, immedi- 
ately upstream. The mass flow rate remains constant at Wj throughout the compres- 
sors and intercoolers. 

For any compressor 


Pc, ex " Pc, in^C 

where the pressure ratio r^ for each compressor is specified by the input variables. 
Either an overall efficiency Vq q ^ polytropic efficiency tiq p can be specified for 
the compressor. With overall efficiency, an ideal exit temperature T^ is found 

by iteration from 


Compressor specific work is then found from the efficiency definition 

ih’ ”(^C.ln’’^C.ex,ld’°) 

”c,o 

and compressor exit temperature from 


(9) 


( 10 ) 
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«(Tt.in-Tc,ex>0) = ^''C 


( 11 ) 


With pol3rtropic efficiency, the compressor exit temperature is found from 

’^C.P 

and compressor specific work from 


( 12 ) 


(13) 


For each shaft, the compressor specific works are summed in order to obtain the 
turbine power required to drive the compressors on that shaft: 


^C,j - 


(14) 


Intercooler . - The inlet conditions to each intercooler are the exit conditions from 
the compressor immediately upstream. The exit temperature from each intercooler is 
specified as program input, and the exit pressure is determined from the input pressure 
recovery rjj^,p 


Pint, ex " Pint, in^iNT 


(15) 


The heat rejected in each intercooler can be fovmd as 

Ah^NT = l^(TiNT, ex’ "^INT, in’ 

Recuperator cold side . - Entering the recuperator cold side is the last (highest 
pressure) compressor exit flow having temperature Tg and pressure Pg. Some of this 
compressor- exit flow is bled off for turbine coolant or to account for a leakage loss such 
as in a rotary recuperator. The mass flow rate entering the recuperator is then 

= -Z*T,c 

The leakage flow w^ and the coolant flow ^ required for each turbine are specified 
by the program input. The recuperator cold-side exit pressure is determined from the 
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input value of pressure recovery r 


R, cold’ 


P3 = P2’’r, cold 

At this point in the analysis it is necessary to enter an iteration loop. The recuper- 
ator cold-side exit temperature Tg cannot be determined imtil the hot-side inlet tem- 
perature Tg is known. The value of Tg depends on the turbine flow rate, which in 
turn depends on the fuel added in the primary burner which in turn depends on the 

temperature Tg. Since the fuel flow is small as compared with the air flow, conver- 
gence is rapid. As a first estimate, 



( 19 ) 


where the primary burner exit temperature T^ is a program input. The remainder of 
the recuperator analysis is presented later in proper sequence. 

Burner . - The burners and turbines are analyzed in flow sequence. Gas entering 
the primary burner comes from the cold side of the recuperator, and gas entering any 
reheat burner comes from the turbine immediately upstream. In the general case, gas 
enters the burner with temperature Tg pressure Pg .j^, and mass flow rate Wg 
having a composition corresponding to a fuel-air ratio of fg For the primary burn- 

'^B,in = PB,in = P3’ "^B,in = ^^ 2 ’ ^B,in = Gas leaves the burner with 
temperature Tl, which is a prc^ram input, and pressure pU which is deter- 

D f GX Jd^ gx 

mined from the input value of pressure recovery rg 


PB,ex = PB,in^B 


( 20 ) 


and with flow 


^B, ex " '^B, in ^f, B 

where Wj g is the fuel flow added to this particular burner, 
ratio is 


( 21 ) 

The associated fuel-air 


^ _’"f,B, ex '^f, B, in ^ ""f, B , , 

^B,ex-:;;^ -'B,in-^^B 

"'air, B '"air, B 


( 22 ) 


where w, „ ;. is the amount of fuel added in previous burners and w„.^ „ is the 
f, B, in ^ air, B 

burner air flow (fraction of original inlet air that exists in the flow in this burner). 

The fuel flow w^ g is found from a burner energy balance, which can be expressed, 
with respect to a reference temperature T^. as 
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"f. BHf(Tr" Tf) + Vf, B ^>'cb + *B, in»(Tr- T'b, in' 'b, in) 


= »'B,exH(Tr.TB,ex'%,ex) <23) 

The heat of combustion Ah^^^ must be that at the reference temperature and is a 
program input aloi^ with the burner efficiency The fuel enthalpy term is a 

polynomial expression defined by equations (C2) and (C4) of appendix C. Fuel tempera- 
ture Tj is either the fuel inlet temperature Tj or the fuel-compressor exit temper- 
ature Tj ^ as obtained from equation (C8). The right-hand side (RHS) of equa- 
tion (23) can be expressed in the form of equation (B9) 


RHS = w 


ex L, ^ ['1 * Tb, ex' «) + ex) 

\ B, 6X / 


(24) 


Substituting for f„ from equation (22), equating the constant air flow in terms of 
inlet and exit conditions such that 


w 


w_. 


B, ex 


w 


B, in 




(25) 


and expressing w^ g as 


yield for equation (24) 


"f,'B = %''air,B 


* ■"’“(’’r- '^B, ex' ») * 'b, lnl(Tr' ex) 

*’*'f,B*(Tr'Tk,ex) <2«) 

In accordance with equation (B9), equation (26) can be expressed as 

RHS = i„H(Te' Tb, ex' 'b, in) " »f. B'(Tr' ex) <”) 

Now, substituting equation (27) back into equation (23), combining terms, and rearrang- 
ing yield 
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( 28 ) 



^B,in 


in’ "^B, ex’ in) 

1, 

Hfl 

(Tr.Tj) 

1 + 

% ^^b - I' 

(Tr> Tb, 

ex) 


Turbine . - Flow into any turbine comes either from a burner or from another tur- 
bine. The exact nature of the turbine calculation depends on which shaft the turbine is 
on. For a multishaft system the drive power for all shafts other than the low-pressure 
shaft (j # is 

PT.j=Pc,j’^P W 

where ^ is obtained from equation (14) and Kp is a factor (program input) to allow 
for output power from the shaft. With the shaft power Pp j known, the power Pp for 
each turbine is obtained from program input. Turbine specific work is then determined 
as 

Pt, 

Ah’ =— !— (30) 

'^Tjin 

and turbine exit temperature from 

^('^T,ex’'^T,in’^T,in) " ^^T 

Either an overall efficiency Vt n ^ polytropic efficiency 77 ^, can be specified 
for each turbine. With overall efficiency, the turbine ideal work is obtained as 


^^T, id “ 


Ah' 




and turbine ideal exit temperature from 

^('^T, ex, id’ "^T, in’ ^T, in) = ^^T, id 
Turbine exit pressure is then 

Pt, ex " Pt, in ~ ^(^T, in’ "^T, ex, id’ ^T, in) 
With polytropic efficiency, turbine exit pressure is 


Pt, ex “ PT,in 


p 


^('^T,in’ '^T,ex’ ^T, in) 


(32) 


(33) 


(34) 


(35) 
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For a single-shaft system and for the low-pressure shaft of a multishaft system, the 
pressure ratio across the shaft is determined from the last turbine exit pressure 



1 

*“r, hot^D 


(36) 


where the recuperator- exit total-to-static pressure ratio P^/Pq, the recuperator hot- 
side pressure recovery rj^ and the diffuser pressure recovery rj^ are program 
inputs. The individual turbine inlet and exit pressures p)p and pip are then ob- 
tained from the input pressure ratio distribution. With polytropic efficiency specified, 
turbine exit temperature Tip specific work Ahip, and power P,p are determined 
from eqxiations (35), (31), and (30), respectively. With overall efficiency specified, 
turbine specific work AhIp is determined from equations (34), (33), and (32) in that 
order. Turbine exit temperature TL and power are then determined from 
equations (31) and (30), respectively. 

For the analysis model, it is assumed that the coolant flow for each turbine bypasses 
that turbine and then mixes with the exit flow from that turbine. The coolant does not 
contribute to either turbine work or loss , To make any better assumption would require 
detailed turbine design information. If such were available, the turbine efficiency could 
be adjusted to accoimt for effects of the coolant flow. It is assumed that the total pres- 
sure pIp after mixing is equal to the turbine exit total pressure pL . The 

mixed total flow rate is equal to 


^T,ex,m - "^T,in ^T, c 


(37) 


and the air flow is 


w„, 


= W„. rr, ._ + Wr 


air, T, ex, m air, T, in T, ck^ 


(38) 


m; 


The fuel-air ratio is then 


w . \ 

f » air, T, in \ 

T,ex,m T,inl 1 

' air, T, ex, m/ 


(39) 


The total temperature tL after mixing is found from a heat balance for the 

A ^ CiAj 1X1 

mixing process; 


^T, ex, m^('^c’ ^T, ex, m’ ^T, ex, m) ''^T, in^(^c’ *^T, ex’ in) 


(40) 
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The coolant temperature T’ is either the last-compressor exit temperature Tg or 
some input value. 

Diffuser, - A diffuser is included in the system to allow for a loss when diffusir^ 
the high-velocity turbine-exit flow to the low velocity desired in the recuperator. The 
diffuser- exit total pressure is determined from the input value of pressure recovery rj^ 

P'e'PsV <“> 


All other variables retain the same values as at the exit of the last turbine. 

Recuperator . - The recuperator cold-side flow Wg was defined by equation (17) 
and the cold-side exit pressure Pg by equation (18). The hot-side flow is the diffuser^ 
exit flow Wg, and the hot-side exit total pressure pl^ is determined by the input value 
of pressure recovery rj^ 

P? " Pe^R, hot 

The cold-side exit total temperature Tg, which was initially estimated using equa- 
tion (19), can now be computed using the recuperator effectiveness definition; 

h(t'2,TJ,0)=.^h(t^,T^,0) (43) 

If the computed value of Tg does not agree with the initial estimate or previously com- 
puted value, the calculation loops back to the primary burner for another iteration. 

After a satisfactory value of Tg is found, the hot-side exit total temperature is obtained 
from a recuperator heat balance: 

WgH(T;,,T^,f6) = W2H(T^,T^,0) (44) 

If there is a recuperator leakage flow, such as in a rotary recuperator, the recu- 
perator hot-side exit throughflow mixes with the leakage to give an exit flow rate of 

Wg=W7 + w^ (45) 

The mixed temperature Tg is then obtained from a heat balance for the mixii^ process; 

T^. fg) = w,h(t^, t;,, f,) (46) 


Overall Performance 


The net shaft output power is 
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( 47 ) 


n 


sh 


^sh,net " (^T,j " ^C,j^ " C 

j-1 

where the fuel-compression specific work ^ is obtained from equation (C9). The 
net plant output power is 


P = 7 ? P 

net 'cv sh, net 


(48) 


where the conversion efficiency reflects a generator, gearbox, or other device. 
The specific fuel consumption is 


3600 w^ ^ ^ 
SFC = 

^net 


(49) 


and the cycle efficiency is 


net 


'ey 


w 


f,tot ^^cb 


(50) 


DESCRIPTION OF INPUT AND OUTPUT 

This section presents a detailed description of the program input, normal output, 
and error messages. The input and corresponding printed output for an example power 
system are included for illustrative purposes. 


Input 

A general description of the program input is given and then followed by an illustra- 
tive example. 

General input . - The data and option indicators are input in data records havir^ the 
NAMELIST name INPUT. The variables and indicators that compose INPUT along with 
descriptions, units, and special remarks are presented in the list to follow. Either SI 
units or U. S, customary units may be used with this program. In this list the symbol- 
ism (X(J), J = 1, N) means that X is a singly subscripted variable having N values 
(X(l) to X(N)) to be entered into the program unless otherwise indicated. Similarly, 
the symbolism ((X(I, J), I = 1,N)I = 1,M) means that X is a doubly subscripted vari- 
able having N xM Values to be entered into the program unless otherwise indicated. 
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The single subscript refers to the shaft number. The first of the double subscripts re- 
fers to the component number on each shaft, and the second subscript refers to the shaft 
number. The order of subscripting is shown in figure 5(a). All subscripted variables 
are dimensioned to allow for a maximum value of 5 for each subscript. 

Values for some of the variables in the input list are internally preset by the pro- 
gram before reading the input. These internally preset values are used by the program 
if alternate values are not specified by the input. Thus, if a preset value is appropriate, 
that particular variable does not have to be specified in the input. These internally pre- 
set values are shown in the input list that follows: 


Variable 

TSO 

PSO 

W 

RIO 

NS HAFT 

(NCOMP(J), 

J=1,NSHAFT) 

((ETAC(I,J), 

I=l,NCOMP(J)) 

J=1,NSHAFT) 

IE TAG 


KPOLY 


RCMIN 


Description 

ambient temperature, K (°R) 

2 

ambient pressure, N/cm (psia) 
ambient air absolute humidity 
inlet total pressure recovery 
number of shafts 

niunber of compressors on each shaft 

compressor efficiency (See lETAC and 
KPOLY. Only ETAC(1, 1) need be 
input if IETAC=0, ) 

compressor efficiency value indicator: 

0 - all ETAC(I, J)-ETAC(1, 1) 

1 - each ETAC(I, J) must be input 

compressor and turbine efficiency 
type indicator: 

0 - overall (isentropic) efficiency 

1 - polytropic efficiency 

minimum value of overall compression 
ratio 


RCDEL increment in overall compression ratio 

RCMAX maximum value of overall compression 

ratio (Calculations made at incre- 
ments of RCDEL for each compres- 
sion ratio starting at RCMIN and 
ending at RCMAX. ) 


Preset value 
518.7 
14.696 
0.0 
1.0 
1 

1 (for J=1 only) 


0 


1 
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Variable 

(RCSHSP(J), 

J=1,NSHAFT) 

((RCCOSP(I, J), 
I=l,NCOMP(J)) 
J=1,NSHAFT) 

((ICOOL(I, J), 
I=1,NC0MP(J)) 
J-1,NSHAFT) 

((RINT(I,J), 

I=1,NC0MP(J)) 

J-1,NSHAFT) 

IRINT 


((TINT(I,J), 

I=l,NCOMP(J)) 

J=1,NSHAFT) 

ITINT 


WLAOWA 

(NTURB(J), 

J=1,NSHAFT) 

((IBURN(I, J), 
I=1,NTURB(J)) 
J=1,NSHAFT) 


Description 

shaft compression ratio factor, fractional 
power of overall compression ratio 
(Sum of RCSHSP(J) must equal 1.0.) 

compressor compression ra;tio factor for 
each shaft, fractional power of shaft 
compression ratio (Sum of RCCOSP(I, J) 
must equal 1 . 0 for each J. ) 

intercooling indicator: 

0 - no intercooling 

1 - intercooling 

intercooler total pressure recovery (Input 
only for those I, J having ICOOL(I, J)=l. 
See IRINT. Only RINT(1, 1) need be 
input if IRINT=0. ) 

intercooler total pressure recovery value 
indicator: 

0 - aU RINT(I, J)=RINT(1, 1) 

1 - each RINT(I, J) must be input 

intercooler exit temperature, K(*^R) (Input 
only for those I, J having ICOOL(I, J)=l. 
See ITINT, Only TINT(1, 1) need be 
input if ITINT=0, ) 

intercooler exit temperature value 
indicator: 

0 - all TINT(I, J)=TINT(1,1) 

1 - each TINT(I, J) must be input 

recuperator leakage flow, fraction of 
inlet flow 

number of turbines on each shaft 

reheating indicator: 

0 - no reheating 

1 - reheating 


Preset value 
1.0 (for J=1 only) 

1.0 (for 1=1, J=1 only) 
0 (for all I, J) 

0 


0 

0.0 

1 (for J=1 only) 

1 (for 1=1, J=1 only) 
0 (for all other I, J) 
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Variable 

((ETABd, J), 

1=1, NTURB(J)) 
J=1,NSHAFT) 

lETAB 

((RBURN(I, J), 
I=1,NTURB(J)) 
J=1,NSHAFT) 

IRBURN 


((TTI(I, J), 

I=1,NTURB(J)) 

J=1,NSHAFT) 

ITTI 


HVF 

TR 

HOC 

ITF 
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Description 

burner efficiency (Input only for those I, J having 
IBURN(I, J)=l. See lETAB. Only ETAB(1, 1) 
need be input if IETAB=0. ) 

burner efficiency value indicator: 

0 - all ETAB(I, J)=ETAB(1, 1) 

1 - each ETAB(I, J) must be input 

burner total pressure recovery (Input only for 
those I, J havii^ IBURN(I, J)=l. See IRBURN. 
Only RBURN(1, 1) need be input if IRBURN=0. ) 

burner total pressure recovery value indicator; 

0 - all RBURN(I, J)=RBURN(1, 1) 

1 - each RBURN(I, J) must be input 

burner exit (turbine inlet) temperature, K (°R) 
(Input only for those I, J having IBURN(I, J)=l. 
See ITTI. Only TTI(1, 1) need be input if 
ITTI=0. ) 

burner exit temperature value indicator: 

0 - all TTI(I,J)=TTI(1, 1) 

1 - each TTI(I, J) must be input 

lower heating value of fuel at temperature TR, 
J/kg (Btu/lb) 

reference temperature for fuel heating value, 

K (°R) 

mass ratio of hydrogen to carbon in fuel 

indicator for fuel temperature entering burner: 

0 - fuel enters burner at temperature TR 

1 - fuel enters burner at temperature TF 
>10 - fuel enters burner at temperature 

calculated from fuel compression work 
(Second digit indicates number of inter- 
cools back to TFIN during fuel com- 
pression. In this case fuel compression 
power is subtracted from gross shaft 
power output. ) 


Preset value 
0 

0 

0 

18 640. 

760. 

0.16786 

0 
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le 


Description 


Preset value 


temperature of fuel entering burner — 

K (°R) (Input only when ITF=1. ) 

constant in fuel heat capacity equa- — 

tion (eq. (Cl)), J/(kg)(K) (Btu/(lb)(°R)) 

(Input only when ITF > 0. ) 

BF constant Bj in fuel heat capacity equation — 

(eq. (Cl)), J/(kg)(K2) (Btu/(lb)(°R^)) 

(Input only when ITF > 0. ) 

CF constant C£ in fuel heat capacity equation — 

(eq. (Cl)), J/(kg)(K^) (Btu/(lb)(°R^)) 

(Input only when ITF > 0. ) 

TFIN temperature of fuel entering fuel compres- — 

sor, K (°R) (Input only when ITF > 10. ) 

PRFIN ratio of pressure of fuel entering fuel com- 1.0 

pressor to ambient pressure (Input only 
when ITF > 10. ) 

ETACF fuel compressor polytropic efficiency — 

(Input only when ITF > 10. ) 

MWF fuel molecular weight (Input only when — 

ITF > 10.) 

((ETAT(I, J), turbine efficiency (See lETAT and KPOLY. — 

1=1, NTURB(J)) Only ETAT(1, 1) need be input if 
J=l, NSHAFT) lETAT = 0. ) 

IE TAT turbine efficiency value indicator: 0 

0 - all ETAT(I, J)=ETAT(1, 1) 

1 - each ETAT(I, J) must be input 

(POWFAC(J), output power factor for all shafts other 1.0 (for all J) 

J=l, NSHAFT-1) than low-pressure shaft (J=NSHAFT), 

ratio of total shaft power to compressor 
power (Not used when NSHAFT=1. ) 

((TSPLIT(I, J), turbine work factor for each shaft (For 1.0 (for 1=1, J=1 only) 

1=1, NTURb(J)) shafts other than low-pressure shaft, 

J=l, NSHAFT) TSPLIT specifies power of each turbine 

as fraction of shaft power; for low- 
pressure shaft (J=NSHAFT), TSPLIT 
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Variable 

Description 

Preset value 


specifies pressure ratio for each turbine as 
fractional power of shaft pressure ratio. Sum 
of TSPLIT(I, J) must equal 1. 0 for each J. ) 


((WCAOWA(I, J), 
I=1,NTURB(J)) 
J=1,NSHAFT) 

turbine coolant flow, fraction of inlet flow 

0.0 (for all I, J) 

ITCOOL 

turbine coolant temperature indicator: 

0 - turbine coolant temperature equals last 

compressor exit temperature 

1 - turbine coolant temperature equals input 

TCOOL 

0 

TCOOL 

turbine coolant temperature, K (°R) (Input only 
when ITCOOL=l. ) 

— 

R65 

turbine exit diffuser total pressure recovery 

1.0 

ER 

recuperator effectiveness 

0.0 

R32 

recuperator cold-side total pressure recovery 

1.0 

R76 

recuperator hot-side total pressure recovery 

1.0 

RSTEX 

cycle exit static -to-total pressure ratio 

1.0 

ETAETA 

shaft power conversion efficiency 

1.0 

TTOL 

temperature tolerance for iterative calcula- 
tions, K (°R) 

0.1 

KOUT 

output indicator: 

0 - output for each cycle point consists of 

overall performance only 

1 - output for each cycle point consists of 

overall performance plus all internal 
temperatures, pressures, and flow rates 

0 

lU 

units indicator: 

1 - SI units 

2 - U. S, customary units 

2 

Input data for NAMELIST input begins with a $ in the second location on a new 
line, immediately followed by the NAMELIST name, which is INPUT for this program, 
immediately followed by one or more blank characters. Any combination of three types 
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of data items may then follow. The data items must be separated by commas. If more 
than one line is needed for the input data, the last item on each line, except the last line, 
must be a number followed by a comma. The first location on each line should always 
be left blank since it is ignored. The end of a group of data items is signaled by a if: 
anywhere except in the first location of a line. The form that data items may take is; 

(1) Variable name = constant, where the variable name may be an array element or 
a simple variable name. Subscripts must be integer constants. 

(2) Array name = set of constants separated by commas where k*constant may be 
used to represent k consecutive values of a constant. The number of constants must be 
equal to or less than the number of elements in the array. This results in the set of con- 
stants being placed in consecutive array elements, starting with the first element of the 
array. 

(3) Subscripted variable = set of constants separated by commas where, again, 
k*constant may be used to represent k consecutive values of a constant. This results 
in the set of constants being placed in consecutive array elements, starting with the 
element designated by the subscripted variable. 

Illustrative example . - The cycle being used as the illustrative example for input 
and output is shown in figure 6. It is a two- shaft recuperated system with four inter- 
coolers and three reheaters. The arrangement of components along with the proper sub- 
scripts for each component are shown in figure 6. Also shown in the figure are all tem- 
peratures, pressures, flows, and component performance parameters required as pro- 
gram input. The fuel is a gas composed of 80 percent (by volume) methane and 20 per- 
cent ethane. The fuel properties required as program input are found in any appropriate 
reference book. It is desired to compute cycle thermodynamic performance for overall 
compression ratios of 5 to 100 in increments of 5. 

The program input for this illustrative example is presented in table I. It is as- 
sumed that the user is familiar with the rules for NAMELIST input, some of which were 
presented previously, and with the order of array storage. A FORTRAN instruction or 
reference manual should be consulted for complete information of this type. Each line 
of the input form shown in table I represents one data card. The output corresponding to 
this sample input is described in the following section. 


Output 

The program output consists of a description of the input and the computed results. 
This section presents normal output. Error message output is described in the next 
section. 

Table II presents the output that corresponds to the input shown in table I for the 
illustrative example. The first part of the output is the description of the input, which 
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is in narrative form. The first line gives the number of shafts and the ambient tempera- 
ture and pressure. The next line presents the inlet and diffuser pressure recoveries and 
the exit static -to-total pressure ratio. The third line gives the recuperator effectiveness 
and pressure recoveries. The next three lines present all the fuel-related input. 

All input parameters associated with the intercoolers, compressors, burners, and 
turbines are then presented in shaft groupings. For each shaft the top line presents the 
shaft number, number of compressors, fraction of overall compression, and number of 
turbines. There is a line for each compressor, giving its fraction of shaft compression 
and its efficiency. If there is an intercooler preceding the compressor, there is a line 
of output preceding the compressor line and giving the intercooler exit temperature and 
pressure recovery. Then there is a line for each turbine on the shaft, giving the turbine 
power or expansion fraction and the turbine efficiency. If there is a burner preceding 
the turbine, there is a line of output preceding the turbine line and givii^ the burner exit 
temperature, efficiency, and pressure recovery. If the turbine is cooled, there is a line 
of output giving the amount of coolant flow. After this information is presented for each 
shaft, there is finally a line of output giving the shaft-power conversion efficiency and 
specifying the nature of the turbomachinery efficiencies. 

After presentation of all the input information, the computed results are then printed. 
Since KOUT = 0 (preset value not altered by input), the computed results that are printed 
include the overall performance only. There is one line of results output for each com- 
pression ratio. The results output includes the compression ratio, specific power output 
(both in kW and hp), specific fuel consumption, cycle efficiency, and fuel-compression 
power. 

If the input had specified that KOUT = 1, the results output would include all internal 
temperatures, pressures, and flow rates in addition to the overall performance. An ex- 
ample of this detailed output is shown in table IE, which gives the results output for a 
compression ratio of 50. There is one Ene of output for each compressor and each tur- 
bine as well as for each side of the recuperator. The Unes of output are in flow se- 
quence. For the compressors and turbines, the shaft and component numbers are given, 
followed by a "NO" or "YES" to specify whether there is an intercooler preceding the 
compressor or a burner preceding the turbine. Then, the next five columns give the in- 
let flow, pressure, and temperature, and the exit pressure and temperature for each 
component. For the turbines the temperature and flow after coolant mixii^ are given 
in the next two columns. This same information is provided for mixing of any leakage 
flow in the recuperator. The last column gives the specific work for each compressor 
and turbine. The overall performance is presented to the right of the table of detailed 
output. 
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Error Messages 


The program contains five output messages indicatii^ the nonexistence of a solution 
satisfying the specified input requirements. These messages are presented in this sec- 
tion, and their causes are discussed. 

(1) COMPRESSOR EXIT TEMP (XXXX. X) GREATER THAN TURBINE INLET TEMP 
(XXXX. X) - This message is caused by the computed exit temperature T^ from the last 
compressor (Cjj) being greater than the primary burner (B^j) exit temperature T^. It 
indicates that the overall compression ratio is too high and/or the burner exit tempera- 
ture is too low. After the message is printed, the program calls for another data set. 

(2) TURBINE -SYSTEM EXIT PRESSURE (XXX. XX) GREATER THAN TURBINE 
EXIT PRESSURE PTT0(I, J)=XXX. XX - This message is caused by the exit pressure for 
any turbine in the system being less than the available pressure Pg at the last turbine 
exit. It indicates that the overall compression ratio is too low for the cycle to be self- 
sustaining in view of the compressor and turbine inefficiencies and the pressure losses 
in the other components. If the compression ratio is not particularly low, one of the in- 
put efficiencies or pressure recoveries could be excessively low. After the message is 
printed, the program increments the overall compression ratio to the next higher value 
and proceeds with the new calculation. 

(3) RECUPERATOR HOT GAS INLET TEMP (XXXX.X) COLDER THAN COLD GAS 
INLET TEMP (XXXX. X) - This message is caused by the last turbine exit temperature 
Tg being lower than the last compressor exit temperature Tg with a recuperator in the 
system. It indicates that the overall compression ratio is too high to provide a temper- 
ature potential for recuperation. After the message is printed, the program calls for 
another data set. 

(4) TURBINE POWER LESS THAN COMPRESSOR POWER ON SHAFT J - This mes- 
sage is caused by the turbine power for the low-pressure shaft (J=NSHAFT) being less 
than the power required to drive the compressors on that shaft. The basic causes for 
this condition are the same as for error message (2), as is the program operation 
sequence. 

(5) OXYGEN USED UP IN BURNER I ON SHAFT J - This message is caused by the 
fuel-air ratio in the flow leaving a burner being greater than the stoichiometric fuel-air 
ratio. It indicates that all the oxygen in the inlet air was used up before all the required 
burning occurred. This can be caused by overall compression ratio being too low or 
burner-exit temperature being too high. After the message is printed, the program in- 
crements the overall compression ratio and proceeds with the new calculation. 
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PROGRAM DESCRIPTION 


This computer program consists of main program MULTI, subroutine FUEL, and a 
set of combustion-gas thermodynamic-property subprograms containii^ the six functions 
CP, H, T2H, S, T2S, and HF and subroutine THERMO. The entire program is written 
in IBM 7090/7094 FORTRAN IV language. In this section, the functions of the main and 
subprograms are described, the program variables are defined, and the program listing 
is presented. 


Main Program MULTI 

Main program MULTI performs all input and output operations, all logic associated 
with the shaft and component arrangements, and all computations except those directly 
using the thermodynamic properties of the fuel and of the combustion gas. 

Program variables. - The variables used in MULTI are defined as follows: 


ABSBF 

absolute value of BF 

ABSCF 

absolute value of CF 

AF 

constant A^ in eq. (Cl) 

BF 

constant B^ in eq. (Cl) 

CF 

constant Cj in eq. (Cl) 

DELHC(I, J) 

compressor specific work 

DELHF 

fuel compression specific work 

DELHIN(I, J) 

intercooler heat removal per pound of flow 

DELHSH(J) 

sum of compressor specific work for shaft J 

DELHT(I, J) 

turbine specific work 

DHID 

turbine ideal specific work 

DHL 

recuperator leakage mixing heat balance term 

DHM 

turbine coolant mixing heat balance term 

DHR 

recuperator heat balance term 

DLHC 

overall sum of compressor specific work 

ER 

recuperator effectiveness 

ETABd, J) 

burner efficiency 

ETACd, J) 

compressor efficiency 
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ETACF 

ETACY 

ETAETA 

ETAT(I, J) 

EXPl 

EXP 2 

FA 

FBOA(I, J) 

FLPOKW 

FOATI(I, J) 

FOATOM(I, J) 

FOA5 

FOA6 

FOA7 

FOA8 

FSTOIC 

FUELPO 

H 

HF 

HFF 

HOC 

HVF 

I 

IBURN(I, J) 

IC 

ICOOL(I, J) 
lETAB 
IE TAG 


fuel compressor efficiency 
cycle efficiency 

shaft power conversion efficiency 
turbine efficiency 

data statement word EXPAN for output use 
data statement word SION for output use 
fuel-air ratio 

fuel-air ratio addition in burner 

fuel compression power 

fuel-air ratio at turbine inlet 

fuel-air ratio after turbine coolant mixing 

fuel-air ratio at station 5 

fuel-air ratio at station 6 

fuel-air ratio at station 7 

fuel-air ratio at station 8 

stoichiometric fuel-air ratio 

fuel compression power 

function defined by eq, (B7), see Combustion- Gas Thermodynamic 
Property Subprograms 

function defined by eq, (BIO), see Combustion-Gas Thermodynamic 
Property Subprograms 

fuel enthalpy relative to reference temperature 
mass ratio of hydrogen to carbon in fuel 
lower heating value of fuel 
dummy index 

reheating indicator, see Input section 

compressor number index 

intercooling indicator, see Input section 

burner efficiency value indicator, see Input section 

compressor efficiency value indicator, see Input section 


lETAT 

IRBURN 

miNT 

IS 

IT 

ITCOOL 

ITF 

ITINT 

ITTI 

lU 

J 

JJ 

JJJ 

KOUT 

KPOLY 

L 

M 

MINUS 

MW 

MWF 

N 

NCOMP(J) 

NS 

NS HAFT 

NTURB(J) 

OFF 

ON 

OVHP 

OVPOKW 


turbine efficiency value indicator, see Input section 
burner total-pressure recovery value indicator, see Input section 
intercooler total-pressure recovery value indicator, see Input section 
shaft number index 
turbine number index 

turbine coolant temperature indicator, see Input section 
burner-inlet fuel temperature indicator, see Input section 
intercooler- exit temperature value indicator, see Input section 
burner-exit temperature value indicator, see Input section 
units indicator, see Input section 
dimensional constant 

turbine number index for last turbine on shaft 
munber of turbines on shaft number NSHAFT 
output indicator, see Input section 

compressor and turbine efficiency indicator, see Input section 
fdummy index 
Iturbine number index 
dummy index 

data statement word for minus sign for output use 

combustion-gas molecular weight function, eq. (B15) 

fuel molecular weight 

index for number of turbines on shaft 

number of compressors on each shaft 

number of shafts minus one 

number of shafts 

number of turbines on each shaft 

data statement word NO for output use 

data statement word YES for output use 

net output power 

net output power 
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PHP 
PKW 
PLUS 
PNC 2 

POWFAC(J) 

POWSH(J) 

POWT(I,J) 

POWTSH 

POWl 

POW2 

PSI 

PRFIN 

PRSH 

PRT(I, J) 

PSC 

PST 

PSO 

PT 

PTCKI, J) 
PTCO(I, J) 
PTTI(I, J) 
PTTO(I, J) 
PTTOM(I, J) 
PTO 
PTl 
PT2 
PT3 
PT5 
PT6 
PT7 


data statement word HP for output use 

data statement word KW for output use 

data statement word for plus sign for output use 

data statement word N/CM2 for ou^Jut use 

shaft output power factor, see Input section 

total power for each shaft except for shaft number NSHAFT 

shaft power for each turbine 

total power for shaft number NSHAFT 

data statement word POW for output use 

data statement word ER for output use 

data statement word PSIA for output use 

ratio of fuel inlet pressure to ambient pressure 

total turbine expansion ratio for shaft nvunber NSHAFT 

expansion ratio for each turbine on shaft number NSHAFT 

value of function $ in eq. (9) or (12) 

value of function $ in eq. (34) or (35) 

ambient pressure 

total pressure 

compressor inlet total pressure 

compressor exit total pressure 

turbine inlet total pressure 

turbine exit total pressure 

turbine exit total pressure after coolant mixii^ 

total pressure at station 0 

total pressure at station 1 

total pressure at station 2 

total pressure at station 3 

total pressure at station 5 

total pressure at station 6 

total pressure at station 7 



PT8 

QF 

R 

RBURN(I, J) 

RBURNT 

RC 

RCCOSP(I, J) 

RCDEL 

RCMAX 

RCMIN 

RCOMP(I, J) 

RCSHFT(J) 

RCSHSP(J) 

RINT(I, J) 

RSTEX 

RIO 

R32 

R65 

R76 

S 

SFC 

SHPOBT 

SIGNBF 

SIGNCF 

TCOOL 

TEST 

TF 

TFIN 

TINT(I, J) 


total pressure at station 8 . 
total heat input to cycle • 
gas constant 

burner total pressure recovery 

product of all RBURN(I, J) for shaft number NSHAFT 
overall compression ratio ■ 

compressor compression ratio factor for each shaft, see Input section 

increment in overall compression ratio' 

maximiun value of overall compression ratio 

minimum value of overalTcompression ratio 

compression ratio for each compressor 

compression ratio for each shaft 

shaft compression ratio factor, see Input section 

intercooler total pressure recovery 

cycle-exit static-to-total pressure ratio 

inlet total pressure recovery 

recuperator cold-side total pressure recovery 

diffuser total pressure recovery 

recuperator hot-side total pressure recovery . 

function defined by eq. (B12), see Combustidn-Gas Thermodynamic 
Property Subprograms 

specific fuel consumption 

net output shaft power 

output word set equal to PLUS or MINUS as appropriate 
output word set equal to PLUS or MINUS as appropriate 
turbine coolant temperature 

value used to test for maximum compression ratio 
temperature of fuel entering burner 
temperature of fuel entering fuel compressor 
intercooler- exit temperature 
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TPK 

TPKl 

TPK2 

TPR 

TPRl 

TPR2 

TR 

TSPLIT(I, J) 

TSO 

TT 

TTCI(I, J) 

TTCO(I, J) 

TTC0ID(I, J) 

TTI(I, J) 

TTOL 

TTTI(I, J) 

TTTO(I, J) 

TTT0ID(I, J) 

TTTOM(I, J) 

TTO 

TTl 

TT2 

TT3 

TT3PRE 

TT5 

TT6 

TT7 

TT8 

TYPEF 

TYPEFP 


data statement word K for output use 

data statement word KELV for output use 

data statement word IN for output use 

data statement word R for output use 

data statement word RANK for output use 

data statement word INE for output use 

reference temperature for fuel heating value 

turbine work factor for each shaft, see Input section 

ambient temperature 

total temperature 

compressor-inlet total temperature 
compressor-exit total temperature 
compressor- exit ideal total temperature 
burner-exit temperature 

temperature tolerance for iterative calculations 

turbine-inlet total temperature 

turbine- exit total temperature 

turbine-exit ideal total temperature 

turbine-exit total temperature after coolant mixing 

total temperature at station 0 

total temperature at station 1 

total temperature at station 2 

total temperature at station 3 

previous value of total temperature at station 3 

total temperature at station 5 

total temperature at station 6 

total temperature at station 7 

total temperature at station 8 

output word set equal to TYPEFP or TYPEFI as appropriate 
data statement word POLY for output use 



TYPEFI 

T2H 

T2S 

W 

WAIRT 

WAIR2 

WAER5 

WAIR6 

WAIR7 

WAIR8 

WARTOM(I, J) 
WBTU 

WCAOWA(I, J) 

WCAWA 

WF(I, J) 

WFTOT 

WGM 

WJLS 

WKG 

WLAOWA 

WLB 

WT 

WTI(I, J) 
WTOM(I, J) 

WO 

W1 

W2 

W3 


data statement word ISEN for output use 

temperature satisfyii^ function H value, see Combustion-Gas Thermo- 
dynamic Property Subprograms 

temperature satisfying function S value, see Combustion- Gas Thermo- 
dynamic Property Subprograms 

ambient air absolute humidity 

airflow rate 

airflow rate at station 2 

airflow rate at station 5 

airflow rate at station 6 

airflow rate at station 7 

airflow rate at station 8 

• ; ' • -I 

airflow rate at turbine exit after coolant mixing 

data statement word BTU for output use 

coolant flow fraction to each turbine 

total coolant flow fraction 

fuel flow rate to each burner ; 

total fuel flow rate 

data statement word GM for output use 
data statement word JLS for output use 
data statement word KG for output use 
recuperator leakage flow fraction 
data statement word LB for output use 
gas flow rate 

gas flow rate at turbine inlet 

gas flow rate at turbine exit after coolant mixing 

gas flow rate at station 0 

gas flow rate at station 1 

gas flow rate at station 2 

gas flow rate at station 3 
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W5 gas flow rate at station 5 

W6 gas flow rate at station 6 

W7 gas flow rate at station 7 

W8 gas flow rate at station 8 

XI output word set equal to POWl or EXPl as appropriate 

X2 output word set equal to POW2 or EXP2 as appropriate 

X3 output word set equal to TPRl or TPKl as appropriate 

X4 output word set equal to TPR2 or TPK2 as appropriate 

X5 output word set equal to TPR or TPK as appropriate 

X6 output word set equal to PPSI or PNC2 as appropriate 

X7 output word set equal to WLB or WKG as appropriate 

X8 output word set equal to PHP or PKW as appropriate 

X9 output word set equal to WLB or WGM as appropriate 

XlO output word set equal to WBTU or WJLS as appropriate 

YN output word set equal to ON or OFF as appropriate 

Program listing . - The FORTRAN listing for main program MULTI is as follows: 


c 

c thermodynamic cycle analysis for MIJLTISHAFT power system with 
C MULTIPLE INTERCOOLS AND REHEATS 

C 

COMMON/CFUEl /AF ,BF,CF tTFIN, TF.ETACF ,I TF,RC,H= F , TR , D ELHF ,MW F , PRF IN 
COMMON FOC,TTOL,W 
REAL J, MW, MWF .MINUS 

DIMENSION WCAOWAI 5,5) , I COOL ( 5 , 5 ) , I BURN ( 5 , 5 ) ,R CC CS R ( 5, 5 ) ,TS PL I T ( 5,5 
1) ,RBURN(5,5), TTCO( 5,5) ,RTC0(5 *5) , TI NT ( 5 , 5 ) , TT C I ( 5 , 5 ) , PTCK 5 , 5 ) , 

2R INTI 5, 5 ) ,DELHIN( 5,5) ,PCQMP (5,5) , TTCOI D ( 5 , 5 ) , OE LHC ( 5, 5 ) ,ET AC( 5 , 5 ) , 
3TTI ( 5, 5) , TTTOM( 5, 5) ,PTTOM (5 ,5) ,WTOM(5,5 ) , WART OM (5 , 5 ) , FOAT 0M( 5, 5 ), 
4PTTK 5, ‘ ) ,TTT I ( 5, 5) ,F30A( 5,5) ,E TAB (5,5) ,F(nATI (5,5) , WF ( 5 ,5 ) , WT I ( 5, 5 
5) ,POWT ( 5,5) ,DELHT(5,5 ) ,ETAT(5,5) ,TTTOIO (5 ,5 ) , PTTD( 5 ,5 ) , PRT ( 5, 5 ) , 
6TTTO( 5,5) ,POWFAC( 5) ,NCONR( 5) ,NTUR3( 51 ,RCSHSP( 5 ) ,FCSHFT ( 5), DELHSHI 5 
7) ,POWSF ( 5 ) 

C 

NAM EL IST/IMPt)T/TS0,PS0,R10,ETAC ,WCA0WA,ER ,R32 ,HVF, T R, ET AR, R INT , 
IRBURM ,TTI , ETAT,FTAETA, W, ICnOL ,I BURN, I TC OOL , I E T AC , I ETAT , lET AB, 

2P0W FAC, N SHAFT, NC0MO,'JTURa,RCSHSP,RCC0SP,TSPLl T, ITTI , I T 1 NT , HOC , R (55, 
3R7iS,T'TCL, TINT,RCMTM,RCntL ,R CMA X , TCOOL , TF , AF , B F , CF , ETAC F, MW F , TF I N, 
AIR BURN, IP IMT, I TF , KO')T , KPQ LY , WLA CWA , PRF IN, PS TEX, ID 
MW( F) = ( 1 .0+F4-W) / ( ,03A52 2<- (F /( l.O+HOC) ) « ( . 2A8D 2*HDC ) i-W / IB.3 16) 

R( F )= 15A5./MW(R ) 

DATA PLlJS,MINI)S/lH+,lH-/ 

DATA PCW1,P0W2,EXP1,E XP?/5H P0w,5HEP , 5HEX PAN, 5HSI ON f 
data TYRFFP, TYPEF I /5H POLV.SH I SEN/ 
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OAT A [JN,nFP/5H YtS,PM NO/ 

OAT A TCP l,TP? 2t TPK 1 ,T PK2 t T PR , TPK /AH HANK , 3HI HP , A L V, 3 H ! 0 , IHR, IHK 

1/ , nGM , RPTU, W.l LS/2HGM, 3HR TH , 3H JL S / 

DATA PPSI,PNC2/6H PSIA , AH N/CH2/ 

DATA WLR, WKG/2HLB , 2HKG / ,PHP ,PKW/2HHP, 2HKV</ 

C 

C INITI AL I/ATIGN 

C 

TSO = 5U.7 
P S0 = 1A .696 
W=0.0 

tTAETA=1.0 
R 10 = 1 . 

DO 72 L=l,5 
DO 71 1,5 

WCAQWA(L,M)=J.O 
ICOOL(L,M ) = 0 
71 inUPN<L,M)= 0 
7 2 PTvJFAC ( LI =1.0 
NSHAFT=1 
NC 3MP f ] )= 1 

\mjPB(l)= 1 
R C S HS P ( 1 I = 1 . 

Rccnsp (1,11= 1 . 

TSPLI T( I, 1)= 1. 

J.RHPN ( 1, 1 )= 1 

ITT1=C 

ITINT=C 

ITCOOL=C 

ItTAC=C 

IETAT=0 

IgTAR=C 

IR INT=C 

IRRUP.N=C 

KnUT=C 

KPOLY= 1 
WLA0WA=0.0 
ITF = 0 
ER=C.O 
R32 = l. 

Hnc = . 16 786 
HVF=18£ AO. 

TR=76C. 

PRF IM= 1 .0 
P65 = l. 

R76=l . 

PSTEX=1. 

IU=2 
.1=778 . 

TT0L=.1 

1 IF( lU .EQ.2) GO TO 73 
TSG=TSC/1.8 
on 76 I=1,NSHAFT 
II=NCnYP(II 
no 7A L=lfII 

7A TINKLt 1)=TINT(L,I 1/1.8 
IJ=NTUPB( I ) 

DO 75 L = lfr.l 

75 TT I (L , I )=TTI ( L, I ) /l.R 
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76 CONTIMLF 

HVP=HVF* loss. P7/. ^5 26 <52 3 7 

TR = 7R/l.a 

TF=TP/l.a 

AF = AF=(= inSS.37/ . AS3S92 3 7<'1. 3 
ar = RF=t=lC55.37/.A53S923 7»!<L.8*L.9 
CF=CF* 10S5.37/.AS3S92 3 7*1. 3#l.8^1.0 
I F( ITF .Gi. .10) TFIM = T'^ INJ/l. 8 
IF( ITGf DL.fcO. 1) TCnnL=TC0nL/1.8 
73 P EAD (5, INPUT ) 

X 3= TP R 1 
X 4= TOP, 2 
X5=TPR 
X6=PPSI 
X7=WI B 
XB=PHP 
X9= WLP 
X10= WRTU 

IF( IU.E0.2) on TIT 7 

X3=TPK1 

X4=TPK2 

X5=TPK 

X6=PNC2 

X7=WKG 

X8=PK W 

X9= WGM 

X10= WJLS 

7 TeST=PCMAX+,l*RCOeL 

FST0IC = .231A«‘ ( l.fHOC) /(2. 6644+ 7. 936 5*H0C » 

IF( ITF.EO.O) TF=TP. 

JJJ=NTUPR(MSHAFT) 

WCAWA=C.O 
PBURMT=1.0 
on 13 P=1,NSHAFT 
L=NCOMP(M) 

IF(L.EO.O) G3 TO 12 
on 11 1=1, L 

IF( icoru I,M) .EO. 0) GO TO 11 

IF( IT INT.IQ.O) TINT (I , M ) = TI N T ( I , I ) 

IF( IP INT.EO.O) RINTd ,V»=PIWT(1 ,1) 

11 IF( lETAC.FO.O » ETACd ,V)=ETAC(1,1) 

12 N=NTUP P( M ) 
no 13 1 = 1,N 

WCAWA = V^CAWA+WCAOWA ( I ,M) 

IF( IPUPN(I,M).EO.O) GO TO 13 

IF( lETAB.EO.O) ETABd , M ) =E TAB ( I , I ) 

IF( IPBUPN.FQ.O) F.BURN ( T ,M ) =PBURN( 1, I ) 

IF( M.EC.NSHAFT) P BUPNT=RBURNT*RRURN(I ,M) 

IF( ITTI.EO.O) TTI (I ,N)=TTI ( 1,1) 

13 IF( IETAT.EQ.O) ETATd ,M)=ETAT( 1,1) 

WRITE INPUT VALUES 

WRITE (6,1000) MSHAFT ,TS0,X3 ,X4 ,PSJ,X6 ,R10 ,P65, 

IRSTEX, f P,P32,P 76, WLAOWA,HVF ,TP , X3,X4,H0G 
1000 FORVATf IHl, I I, 54H-SHA FT POWhP SYSTEM OPtRATING AT AMBIENT CONDITIO 
INS nF,F8.2,S-< neo ,A4,A3,4H AND,F8.3,A6/ 3 1H INLET TOTAL PRESSU 

2RE PEGOVFF Y=, F 5.3, 54-1 , TURBINE EXHAUST DIF FUSER TOTAL PRESSURE 
3REC0VFRY=,P5.3, 29H, EXIT STAT/TCT PRESS FATIO=,R5.3, 



3 /27H RFCtJPIiPATnR tPFSCT I Vc'gPSS= , F“5.3 , 36H, rOLO S I n;- 

-VTOTAL PPFSSU7E PrCH Vc R V= , F 5. 3 , 3 5H, HOT S I Ht TOTAL PRiiSS'JPF RFCOVrR 
5V = , F5.2t IRH, bypass FL0W=,F6.4 

5/20H FLEL HtATINS \ML Ut = , F 9. 0, 3H AT,F6.0,5H DEG ,A4,A3, 12H, H 

6/C P4T10=,F7.5) 

IF( ITF.GE .10) GO TO 42 

41 WRITE (6,1010) TF,X3,X4 

1010 FORMAT ( 1H + , 65X, 23H, FUEL £^iTETS BURNER AT,F6.0,5H OEG ,A4, A3/) 

IF( ITF .EO.O) GO TO 43 

42 SIGNBF=pLLS 
SIGNCF=MINUS 

IF( RF.LT.O.O) SIGNBF=MINUS 
IF(CF.CE.O.O) SIGNCF=PLUS 
ABSBF = ABS(BF ) 

ABSCF=APS(CF) 

WRITE (6,1020) AF,SIGNBF,ABSBF,SIGiMCF,A8SCF 

1020 P0RMAT(13H FUEL GAS S P = ,F 9. 4, 1 X , A1 , 1 PEI 1 . 4 , 5H * T ,A1,E11.4, 7H * 
1T*<'2) 

I F( ITF .LT ,10) GO TO 43 

WRITF (6,1021) MWF.TF IN,PRF IN.cTACF 

1021 FORMAT! 14H 

1 FUEL MOL WGT=,OPF7.3, 17H FUEL INLET TEMP= ,F 5 , 0, 25H, FUEL INLET C 
20MPRf;SSI0N'=,F 5, 1, 24H, Fyt L COMPRESSION £FF. = ,F9.3) 

43 on 50 I=1,NSHAFT 
IS=NSFAFT + 1-I 
L=NCOMP( IS ) 

M=NTURB( IS ) 

WRITE (6,1030 ) IS,L,RCSHSP(IS) ,M 

1030 FORMAT( 7H0SHAFT , I 1 , 5H HAS ,I1,24H COMPRESSORS THAT SUPPLY , F7 .4, 26 
IH OF TOTAL COMPRESSION AND , 1 1 , 9H TURBINES) 

IF( IS.LT.NSHAFT) WRITE (6,1040) POWFAC(IS) 

1040 FORMAT( 1H*,81X, 12HTHAT OE L I VE R , F7. 4 ,2 5H * COMPRESSOR SHAFT POWER) 
IF(L.EO.O) GO TO 46 
DO 45 LL=1,L 
IC=L<-1-LL 
WRITE (6,1050) 

1050 FORMAT( IH ) 

I F( ICOOL( IC, IS ).NE.O) WRITE (6,1060) IC,TINT( IC, IS) ,X3,X4,RINT( IC, IS) 

1060 FORMAT ( 6X,24HFL0W ENTERING COMPRE SSOR, I 2 , 13H IS COOLED T0,F6.0,5H 
lOEG ,A4,A3,30H WITH TOTAL PRESSURE RECD VE RY= , F5 , 3) 

45 WRITE (6,1070) IC , RCC QSP ( I C ,I S) ,ETA C ( I C , I S ) 

1070 F0RMAT(6X,10HCnMPp£SSnR,I2,9H PROVI DE S , F7 . 4, 3 8H OF SHAFT COMPRESSI 
ION WITH EFFICIEMCY=,F5.3) 

46 DO 50 n=l,M 
WRITE (6,1050) 

I F( 1BURM( IT, IS). NE.O) WRITE (6,1080)1 T,TTI I IT ,I S) ,X3 ,X4,ETAB( IT, IS), 
1RRURN( IT, IS) 

lOBO F0RMAT(6X,21HFL0W ENTERING TUP8 I NE , 1 2 , 1 3H IS HE ATE D T 0 , F7. 1 , 5H DEG 
1 ,A4,A3,29H WITH COMBUSTION EFFICIENCY =,F5.3,29H AND TOTAL PRESSU 
2RE RECCVEPY = , F5.3) 

I F( IS-NSHAFT) 47,48,48 

47 X i=pnw 1 
X2 = POW2 
GO TO A9 

48 X1 = EXP 1 
X 2= EXP 2 

49 WRITE (6, 1090) IT , T SP L I T ( I T , I S) , XI , X2 ,£ TAT ( IT , I S ) 

1090 F0RMAT(6X,7HTURBINE,I 2,9H PROVI DE S, F7.4 ,10H OF SHAFT ,2A5,16HWITH 
lEFF IC IENCY=,P 5. 3) 
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IF(WCACW/i( IT, ISl.EO.O.O) GO TO 50 
WPI'^e (6,1100) IT,UCA0WA( I r, I Si 

1100 FORMAT ( 6X,7HT(IRBIME ,1 2, 15H IS COOLLO WI TH ,F7. 4, 24H GF CCHPRESSOR F 
IX IT flow ) 

I F( ITCrOl .NE.O ) WRITE (6,1110) TCOOL,X3,X4 
mo FORMAT( lH+,6lX,9HCnOLen T0,F7.1,5H OtiG ,A4,A3) 

50 CONTINLF 

TYPEF=TYPFFP 

I F( KPOLY.FQ.O) TYPEF=TYPfcFI 
WRITE (6,1120) FTA£TA,TYPcF 

1120 FORMAT! 20H00DTPUT POWER E OUAL S, F7 . 4 , 1 5H OF SHAFT POWER, 20X , 31HTURR 
lOMACHINERY EFFICIENCIES ARE ,A5, 6HTP GPIC ) 

IF( KOdT.FO.O) WRI TE (6,2000) X7,X7,X8 
2000 FORMATiani PRESS. , 6X, 3HME T , 8X , 3HNE T , 8 X , 4HFUE1. , 7X , 5 HCYC LE,6 X , 4HFUEL 
1, 18X, 18HPOWEP VALUES REFER /7H RATI 0 , 6 X , 5HP0 WE R , 6X , 5HP0W6 R, 5X , 8HC 
20NSUMP . ,6X,3HEFF,6X ,5HP0WER ,18X,7HT0 ONE ,A2,8H PEP S EC/ 14X , 2HKW , 9 
3X,2HHP,7X,A2, 4H/HR/,A2,17X,2HKW,19X,17H0F INLET AIR FLOW) 

IF( lU .EO. 2 ) GO TO 16 
TS0=TSC*1.8 
17 no 36 I=1,NSHAFT 
I I=NCOMP( I ) 

00 21 L=l, II 

21 TINT( L, I )=TINT(L, I )*1. 8 
I J = NTURB( I ) 

00 37 L=1,IJ 

37 TTKL, 1 )=TTl(L,n*1.8 
36 CONTINLE 

HVF=HVF/1055.a7*.45359237 
TR = TR* 1.8 
TF=TF* 1.8 

A F=AF/ 1055. 87*. 453592 3 7/1 .8 
BF=BF / 1055. 87*. 453 5923 7/ 1.8 /I. 8 
CF= CF/1055.87*.45359237/ 1. 8/1. 8/1. 8 
I F( ITF .GE .10) TFIN= TFIN*1.8 
IF( ITCnOL.FO.l) TcnOL = TCnOL*1.8 

INLET CONDITIONS 

16 TTO= TSC 
PTO=P SO 
W0= 1, +W 

INL ET 

5 TT1=TTC 
PT1=PT C*R 10 
W1 = W0 
RC = RCM IN 

6 CONTINLE 
OL HC=C .0 
CALL FUEL 

START INTERCOOLFR-COMPRESSOR LOOP 

DO 150 I=1,NSHAFT 
IS= NSFAFT + 1- I 
RCSFFT(IS)= RC**PCSHSP( IS) 

DEL FSH( IS )=0.0 
IF( I- 1 ) 2,2,3 
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2 TT = TT1 
PT=PT1 
GO TO A 

3 TT = TTCO( I, IS+ 1 ) 

PT= PTCQ( 1, IS + l ) 

A L=NCnWF( IS ) 

IF(L.fcG.O) GO TO 
on 14C M = 1,L 

IC = L4-l-M 

IF(M.EC.l) GO TO 131 
TT= TTCO( IC+1, I S) 

PT=PTCC( IC+1. IS) 

INTEPCGOLFP 

131 IF( ICnCK IC, IS ).EQ.O) GO TO 132 
TTC I( IC, IS )= TINT! IC, I S) 

PTCI(IC,IS)= PT*R INTdC ,I S) 
nELHIN(IC,IS)= H(TTCI ( I C , I S ) , TT , 0. 0 ) 

GO TP 13 3 

132 TTC I( IC, IS)= TT 
PTC I( IC, IS )= PT 
OELHIN ( IC , IS)=C.O 

COMPRESSOR 

133 RCOMP( IC, IS) = RCSHFT{IS)**RCCOSP(IC, IS) 

PTCOI IC, IS)= PTC I (IC, IS)*RCOMP(IC,I S) 

psc=R ( c.o j/.if'ALOGiRcnMPn c,i sn 

IF(KPOLY.EO.l) GO TO 13A 

TTCOIOI IC, IS) = T2S(TTCI ( IC ,I S) ,0.0,PSC) 

OELHCI IC, IS) = H( TTCI (IC ,IS) ,TTCOID( IC , I S ) , 0.0) / ETAC ( IC , I S ) 
TTC0( IC, IS)=T2H(TTCI ( IC,I S) ,0.0,0ELHC( IC, IS)) 

GO TO lAO 

134 PSC=P SC/ETAC( IC,IS) 

TTCO( IC , IS)=T2S(TTCI (IC ,I S) ,0.0,PSC) 

0ELHC( IC, IS) = H(TTCI(IC ,rS) ,TTCO(IC,IS) ,0.0) 

140 r)ELHSH(IS)= OELHSH(IS)+DELHC(IC,IS) 

GO TO 150 

141 TTC0( 1, IS )=TT 
PTCt’( 1 , I S )=PT 

150 DLHC=DLFC«-OELHSH( IS) 

TT2=TTCO( 1,1) 

IF( ITCOCL.EQ.OI TC00L=TT2 
IF( TT2.GT.TTI ( 1, I ) ) GO tq iqO 
PT2=P TCO( 1, 1) 

BYPASS AIR 

W2=W1*( l.-WCAWA-WLAOWA) 

WAIR2 = V,2/W1 

recuperator - initial ESTIMATI MADL FQR TT3 

PT3=R 32*PT2 
W3 = W2 

TT3=TT2+EP*(TT4-TT2) /2. 

P T 5 =p S C /F 7f> /R 6 5 /P. STE X 
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START P.IJFNtR-TURB INe LOOP 

161 WFTPT=C.O 
PDkJTSH= 0.0 
DO 17C IS=1,MSHAFT 

PnWSH( IS )= Wl*DELHSH( I S)*PQWFAC (IS) 

IF( IS- I I 171, 171, 17? 

171 TT=TT3 
PT=PT3 
WT = W3 

WAIRT = KAIP;? 

FA=0.0 

r,n TO 173 

172 JJ=NTURB( JS-l) 

TT= TTTOV( JJ, IS-1 ) 

PT=PTTOM( JJ,IS-l) 

WT=WTC^ (JJ, IS- 1 ) 

WAIRT = WAPTaM( JJ, IS-1) 

FA=FOATrM( JJ, IS-l ) 

173 L=NTURB( IS) 

IF( IS .EC) .NSHAFT) PR SH =PT/ PT5* RBURMT 

DO 180 n=l,L 

IF( IT .EQ. 1) GO TO 174 

TT = TTTCM( IT-1, IS) 

PT=PTTC^'r IT-1, IS) 

WT = WTn>'( IT-1, IS) 

WAIRT = RAPTOM( IT-1, IS) 

FA= FOATOM( TT-l,IS) 

174 IF( I6URN( IT, IS ) .£0.01 GO TO 175 

BURNER 

"TT I( IT, IS 1= PT^RPURN (IT, IS) 

TTTK n,IS)= TTK IT, I S) 

F30A( IT,1S) = ( 1.+FA+W)*H(TT,TTTI (IT,IS) ,FA) /(cTABdT ,IS) *HVF-HF(TP, 
ITTT I( IT, IS) )+HFF ) 

GO TO 176 

175 PTTK IT, I S)=PT 
TTTK IT, IS)=TT 
FBOA{ IT, IS )= 0.0 

176 FOATI ( IT, IS)= FA + FBOA ( I T , I S ) 

I F( FOATK IT, I S) .GT.FSTOIC ) GO TO 500 
WF( IT , IS )=FB3A (IT,IS)*WAIRT 
WTI ( IT, I S ) = WT + WF( IT,I S) 

WFTOT = VsFTOT+WF ( IT, IS) 

TUP BINE 

I F( IS. EC. NSHAFT) GO TO 178 
P0WT(1T,IS)= POWSH(IS)*TSPLIT(IT,IS) 

DEL HT( IT, IS)=PO)VT(IT, I S ) / WT I ( I T, I 5) 

TTTO( I TTTs)=T 2H(TTTI ( I T,I 5) , FOATI ( I T, IS) ,-DEL HT ( 1 T , IS ) ) 
IFIKPOLY.FO.l ) GD TO 177 
DHin=DELHT( IT, IS) /ETAT( IT ,I S) 

TTTOIDI IT, IS)= T2H(TTTI (I T,I S) ,FOATI ( IT,IS) ,-OHin) 

PTTO( IT, I S)=pTTI( IT ,I S ) *E XP ( J /R ( FOA TI ( I T , I S ) ) *S (TT T I ( I T , IS ) , 
1TTT0ID( IT, IS) , FOATIUT, IS ) ) ) 

GO TO 182 

177 PTTO(IT,IS)=PTTI(IT,IS)>f=E XP ( J/R ( FOATI ( I T, I S ) ) / ET AT ( IT , IS ) << 
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1S(TTTI( IT,IS),TTTn(IT,IS> ,f=OATI (1T,IS») ) 

GO TO 1 fl 2 

178 PPT(IT,IS)= PRSH#*TSPLIT( IT,I S) 

PST=R ( FGATK IT, IS ) )/J*ALOG( 1. /PRT( I T, I S) ) 

PTTO(IT,IS)= PTTI (IT, 1 S)/oRTt I T,IS) 

IFfKPOLY.EQ.l ) GO TO 179 

TTTOITM IT, IS ) = T2S(TTTI nT,ISI,FOATI(IT,IS),PST» 

DELHT( IT, IS)=LTAT( IT, 1 S)*H(TTTQID( I T, I S) , TTTI ( IT,I S ), FOATI ( IT, IS) I 
TTTO( IT,IS)=T2H(TTTI( IT,I S) ,FOATniT,IS) ,-DtL HT ( IT , IS ») 

GO TO 181 

179 PST=PST«ETAT(IT,IS) 

TTTO( n, IS )=T2S(TTTI ( IT,I S) ,FUAT1 ( I T, IS ) ,PST) 

D£LHT( IT, IS) = H(TTTn(IT,IS) ,TTTI ( I T , I S ) , FO AT I ( I T , I S ) ) 

181 PnWT(IT,lS)= WTl (!T,I S)*DELHT{1T,IS) 

POWTSH=PnWTSH*POWT( IT,I S) 

COOLANT MIXING 

182 PTTOMl IT, IS)= PTTO(IT,IS) 

IF(PTTC(IT,IS).LT.PT5#.9999) GO TO 200 
WTnM(IT,IS)= WTI(IT,I S)^-W1*WCA0WA(IT,IS) 

WARTOMf IT, IS)= WAIRTt-WCAOWAdT.IS) 

FOATOMI IT,IS)= FOATK IT,I S)*WAIRT/NART0M(1T,I S) 

I F( WCAOVsA ( IT, IS I.EQ.O .0) GO TO 183 

DHM=WT1( IT, IS ) /WT0M( I T , I S ) *H ( TCOUL , TTTO ( I T , 1 5 ) , FOAT I ( IT, IS ) ) 

TTTOM( IT,IS)=T2H(TCnOL,FQATOM(IT,IS),DHM) 

GO TO 180 

183 TTTOM( IT, IS)= TTTO(IT,IS) 

180 CONTINUE 
17C CONTINUE 

IF(POWTSH.LT.Wl*QfcLHSH(NSHAFTn GO TO 400 
TT5=TTT0M( JJJ ,NSHAFT) 

PT5= PTTGMf JJ J,NSHAFT ) 

W5*WT0M JJJ,NSHAFT) 

F0A5=F0ATnM( J JJ,NSHAF T) 

WAIF5 = V,ART0M( JJJ,NSHAFT) 

RECOVERY n IFF USER 

TT6=TT5 
PT6=PT5*R65 
F0A6=F0A5 
W6 = W5 

WAIR6=WAIR5 

RECUPERATOR 

IF( ER. .EQ.O.O ) GO TO 24 
IF( TT6.lt .TT2) GO TO 300 
DHP.=ER«H( TT2, TT6,0.0) 

TT3PRE^TT3 
TT3=T2K TT2,0.0,0HR ) 

1F( ABS(TT3-TT3PPE I.GT.TTOU GO TO 161 
TT7=T2F( TT6,FOA6,-W2/W6*OHR ) 

GO TO 25 
24 TT7=TT6 
2 5 W7=W6 

FOA7=FOA6 

WAIR7=\%AIP6 
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PT 7=PT f 76 
C 

C LEAKAGb MXIMG 

C 

PT8=PT7 

W8=W7-*-Wl*WLA0WA 
WAIPR = WAIP 7 + WLAOWA 
FnA8=F0A7«WAIR7/WAIP.8 
IF( WLACVvA .EO.O.C) GO TH 27 
DHL=H( 7T2,TT7,FnA 7) 

TT8=T2F( TT2,FnA8,W7/W8*nHU 
GO TO 31 
27 TTd=TT7 

PERFORMANCc PARAMETERS 


31 SHPORT = C.O 
IF(NSHAFT.EQ.I) go to 33 
NS=NSHAFT- 1 

nn 32 1=1, MS 

32 SHPORT = SHPOBT+PnWSH (I ) # ( 1 .-1 . /POWF AC ( I ) ) 

33 SHPOBT=SFPOBT>PnWTSH-Wl*QELHSH( NSHAFT) 

IF( ITF.LT.IO) GO TO 34 

FUELPn=ViFTOT*DELHF 

SHPOflT=SFPORT-FUELPO 

FLPOK W=FUFLPO«= 1.0 542 
GO TO 35 

34 FLPOKU=C.O 

35 0VPOKW=SHPOBT*1.0542*ETAETA 
0VHP=SFP0PT*1.414 5<‘eTA£TA 

S FC =WF TOT* 3600. /O VHP 
OF=WFTOT=»HVF 
ETACY=$FPn8T/QF*ETAETA ' 

TF( lU .EQ.2J GO TO 59 
OVPPKW= nVPOK W/. 45359237 
OVHP= OVHP/. 45359237 
SFC= WFTOT*3600./OVP0KW 
PLPOKW= PtPOKV^/. 45359237 

WRITE OUTPUT 

59 IF(KOUT.EO.I) GO TO 190 

WRITE (6,2010) RC ,0'/P0KW,0VHP,SFC,ETACY,FLP0KW 
2010 FORMAT! IX, F 7.2, 2F 11.2, 3F1 1.4) 

GO TO £00 
190 continue 

IF( PC .EQ.FCMIM ) WPITP (6,3000) 

3000 FORMAT ( IHl ) 

WR ITE ( 6, 3010) OVPOKW, X 7 ,RC , XI 0 , OVHP , X7 , X6 , X5 , X6 ,X5 , X5 , X 9, S PC, X 7 , X 8 
3010 FORMAT ( 114H0PRESSURE RATIO COOL INLET INLET INLET E 
IX IT EXIT MIX MIX WORK NET P0WEP = ,F7.2, 

24H KW/,A2,4H A I R /F 1 1. 2 , 9X , 68HOK PLOW PRESS TEMP PRFS 

3S TEMP TEMP FLOW ,A3,4H Pt R , 1 3X , 1 H= , F7. 2 , 4H HP/,A2,4 

4H AIR /19X,4H0URN,11X, A6,4X,4HDtG , A 1 , 3X , A6, 4X , 4Hf)E G , A 1 , 4X , 4H0EG , 
5A1, 12X,A2,5H FLOW , 9X, lOHF UE L C JNS=, F7.4 , 1 X, A2 ,4H/HP/ , A2 >. 

K=1 

no 310 I=l,NSHAFT 
I S = NSHAFT 4 l-I 
L=NC0MP(IS) 
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IF(L.EC.O) n:i TO 310 
nn 30 fl v = i,L 

IC=L+1-M 

IF( ID .tc.2) on TO lOl 
TTCI(IC,IS)= TTFI (IC, I '■-) /l.P 
TTCn(ICtIS)= TTClMir, IS) /l.f? 

OELHCI IC, IS)= DFLHC (I C , I S ) * I 055. 87 /453 . 5923 7 
191 YN= OFF 

IF< ICnOU IC,TS).N[ .0) VN=HN 

WO ITF (6,3020) IS ,IC, VN.Wl , PTC I ( I C , I S ) , TTC I ( I C , IS ) , PT CO ( I C , I S ) , 
1TTCP( IC, IS),OFI.HC ( IC, IS) 

3020 FnPf^4T(6H SH6FT,I2,7H - C OMP , I 2 , A5 , F9 . 4 , 4F9 . 1 ,2 1 Y , F 7. 2 ) 

if(k.e:c.i) write ( 6 , 3021 ) t tacy 

3021 FORMAT( 1H + , 103X, lOHCYCLt EFF=,F7.4) 

IF(K.EG.2) WRI T£ ( 6,3022) FLP0KW,X7 

3022 F0PMAT( 1H + , 103X, lOHFUEL PnwR = , F 7. 4 , 4H KW/,A2,4H AIR) 

K=K+1 

308 CnNTINLE 

310 CONTINLF 

IF( IU .EO.2) on TO 311 
TT2= TT2/1.R 
TT3= TT3/1.8 

311 WR ITE ( 6,3030) W 2, P T2, TT2 , P T 3 , TT3 

3030 FORMAT! 12H P FC UPE RATOP , 1 OX , F9.4 ,4F9. 1 ) 

IF(K.eC.2) WRITE (6,3022) FLPGKW,X7 
on 32 0 IS=1,NSHAFT 
L=NTURP(IS) 
on 320 n = l,L 
IF( IU .EQ.2 ) CO TO 312 
TTTI(1T,IS)= TTTI (IT,I S)/1.8 
TTTO(IT,IS)= TTT0(IT,IS)/1.8 
TTTGM( IT, TS)= TTT0M(IT,IS)/1.8 
OEI HT ( IT, IS )= DELHT(IT,IS)*1055. 87/453. 59237 
312VN=0FF 

IF( IBUPN( IT, IS).r!E.C) YN = ON 

320 WRITE (6,3040) I S , I T , YN , W TI ( I T , I S) , PTT I ( I T , IS ) , TTT I ( I T , IS ) , PTT 0(1 T 
I, IS),TTTO( IT, IS),TTTOM (IT,I S) ,WTr.M(IT,IS) ,D6LHT (IT , IS ) 

3040 F0RMAT(6H SHAFT, 12, 7H - T UR B , I 2 , A5 , F9 .4 , 5 F9 . 1 , F9 .4 , FIO . 2 ) 

I F( IU .EO.2) CO TO 313 
TT6 = TT6/1.3 
TT7= TT7/1.8 
TT8= TTE/1.8 

31 3 WR TTE ( 6, 3050 ) W6 ,PT6 , TT6 ,P T7 , TTT , TT8 , W8 
3050 F0RMAT(12H RECUPERATOR, 10X,F9. 4, 5F9.1,F9. 4//) 

CO TO 60C 

100 IF(IU.EO.l) TT2= TT2/1.3 

IF(IU.EO.l) TTI ( 1, 1) = TTI ( 1,1) /1.8 
WR ITE ( 6, 50001 PC , TT2 , TTI ( 1 , 1 ) 

5000 FORMAT! IX, FT. 2, 4X, 

1 22HCCMPRESSUR EXIT TEMP (,F6.l,35H) CREATEP THAN TURBIN 

IE INLET TEMP ( ,F6.1,1H) ) 

IF(IU.EO.l) TTI(1,1)= TTI ( l,l)*1.8 
CO TO 1 

200 WRITE (6,6000) RC , RT5 , 1 T , I S ,PTTO( I T , I S ) 

6000 FORMAT! 1X,F7.2,4X,30HTURBINE"SYSTGM EXIT PRESSURE (,F6.2,42H) CREA 
ITER than TURBINE EXIT PRESSURE PTTO ( , 1 1 , IH , , I 1 , 2H ) = ,f 6 . 2) 

CO TO 600 

300 IF(IU.EO.l) TT6= TT6/1.3 
IF( lU.EO.l) TT2= TT2/1.8 
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WPTTE (6,7000) PC,TT6,TT2 
7000 format ( 1X,F7. 2, , 

I 32HP Ff UPtP ATOR HOT OAS INLLT TFMP (,‘=6.l,35H) CriOER TH 

IAN COLD GAS INLET TEMP (,F6.1,1H)) 

GO TO 1 

AOO i^RITE (6,8000) pr,MSHAFT 

8000 FORMAT! IX, F7.2,AX, AOHTURBINF POWER LESS THAN COM"R'=SSOP POKER ON S 
IHAFT, 12 ) 

GO TO 6C0 

500 WRITE (6,9000) RC,IT,IS 

9000 FORMAT! lX,F7.2,4X,24HnXYGFN USED UP IN 8UPNER,I2,9H ON SHAFT, I2» 
600 CONTINUE 

PC = RC+PCC'F;L 

IF( RC-TEST ) 6,1,1 

ENO 


Subroutine FUEL 

Subroutine FUEL performs the computations involving the fuel thermodynamic prop- 
erties and the fuel compression. 

Program variables . - Certain of the variables transfer between main program 
MULTI and subroutine FUEL by means of labeled common block /CFUEL/. These vari- 
ables, which were defined in the MULTI variable list, are AF, BF, CF, DELHF, 
ETACF, HFF, ITF, MWF, PRFIN, RC, TF, TFIN, and TR. The remaining variables 
in subroutine FUEL are defined as follows: 

CP arithmetic statement function for heat capacity of fuel, eq. (Cl) 

DFSCF derivative of FSCF with respect to TF 

FSCF difference between right- and left-hand slides of eq. (C8) 

H arithmetic statement function for fuel enthalpy difference, eq. (C2) 

NCF number of fuel compressors with intervening intercoolers 

PRC pressure ratio across one fuel compressor 

PSCF value of fimction in eq. (C8) 

S arithmetic statement function defined by eq. (C3) 

T temperature 

TFl previous value of TF 

Program listing. - The FORTRAN listing for subroutine FUEL is as follows: 
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SU'^PnUTIMfc PUEil 
REAL VXF 

FUK /AF,BF,CP,TFIN,7f-,cTACF ,1 TF , R C , HF F , tp , r)cLH<= ,«WF, PRF IN 
F p ( T ) = A F + B F <' T # 2 <• C F T =!« =(= 3 

T2-T l»+RF/2.*( T2**2-T1 *’«'2 ) 4-CF/3. « ( T2* >!' 3-T I «*3 I 
S( T 1, T2 )=AF<-aLOG( T2/T1 UBF* ( T2- T1 ) CF /2 . * ( T2i‘ *2-T 1 * *2 ) 

I F( ITF .LT .10) GO TP 2 

MCF =MCn( ITF, 10)+1 

PRC=( FC/PF FIM )**( 1,/FLQAT(NCF) ) 

I F( RC /PRFIN.L T . 1. 0) PRC=1.0 

PSCF= 15A5./'1WF/778. /E T ACP *A LOG ( PRC ) 

TF=TF IN*tXP(P SCF/CP(TFIN) ) 

1 FSCF=P5CF-S(TFIN,TF) 
nFSCF=-CP< TF) /TF 
TFl=TF 

TF=TF-FSCP/DFSCF 

IF( A3S(TF-TF1 ) .GT.. 1) GO TO 1 

DEL HF = FLGAT(NCF )>f< H ( TF I N , TF ) 

2 HFF=H( TR,TF I 
RFTURM 

END 


Combustion-Gas Thermodynamic -Property Subprograms 


All calculations directly involving the combustion-gas thermodynamic properties 
are performed by subroutine THERMO. In order to use the values computed in THERMO 
directly in arithmetic statements of the callii^ program, there are six function subpro- 
grams used in conjunction with THERMO. Each function subprogram calls on 
THERMO to do one particular computation, and the results of this computation is then 
set equal to the function. The six types of computation done by THERMO and the associ- 
ated fxmctions, with argviments, are as follows: 


Function 

CP(T,F) 

H(T1,T2,F) 

T2H(T1,F,X) 

S(T1,T2,F) 

T2S(T1,F,X) 

HF(T1,T2) 


Computation 

heat capacity c from eq. (B5) 

enthalpy difference Ah from eq. (B7) 

final temperature Tg by iteration from eq. (B7) 

entropy function A.(p from eq. (B12) 

final temperature Tg by iteration from eq. (B12) 

fimction I(Tj, Tg) from eq. (BIO) 


The fvmction arguments, which are known values for each computation, are defined 
as follows; 
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F fuel-air ratio 
T temperature 
Tl initial temperature 
T2 final temperature 
X function value of Ah or Acp 

The call on subroutine THERMO is CALL THERMO(TIN, TOUT, FOA, FUNC, IND) 
where the ai^xunents are defined as follows; 

FOA fuel-air ratio 

FUNC function value 

IND computation-type indicator: 

2 - Ah 

3 - Tg from Ah 

4 - Acp 

5 - T 2 from Acp 

6 - I(T^, Tg) 

TIN initial temperature 

TOUT final temperature 

This set of subprograms can be used with any program requiring combustion-gas 
thermodynamic properties . 

Program variables . - The only variables used in the function subprograms are the 
function and its arguments, all of which have been defined. Certain of the variables 
used in subroutine THERMO are transmitted from main program MULTI through un- 
labeled common. These are HOC, TTOL, and W, which were defined in the MULTI 
variable list. The remaining variables in THERMO, with the exception of the functions 
and ai^iunents defined previously, are defined as follows: 


A 

constant 

Aair in eq. (B5) 

AC02 

constant 

A^q in eq. (B5) 
2 

AH20 

constant 

Ajj Q in eq. (B5) 
2 

A02 

constant 

in eq. (B5) 
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A1 

constant 

Kj in eq. (B5) 

A2 

constant 

Kg in eq. (B5) 

A3 

constant 

Kg in eq. (B5) 

B 

constant 

®air 

BC02 

constant 

B^q in eq. (B5) 
2 

BH20 

constant 

O - 

2 

B02 

constant 

Bp. in eq. (B5) 

C 

constant 

^air ^ 

CC02 

constant 

C^Q in eq. (B5) 
2 

CH20 

constant 

Cr o 
2 

C02 

constant 

Cp. in eq. (B5) 
^2 

D 

constant 

°air 

DC02 

constant 

D^q in eq. (B5) 
2 

DFCN 

derivative of function FCN with 

DH20 

constant 

Djj Q in eq. (B5) 
2 

D02 

constant 

Dp, in eq. (B5) 

E 

constant 

Eair in eq. (B5) 

EC02 

constant. 

E^q in eq. (B5) 

EH20 

constat 

% O 
2 

E02 

■ constant 

Ep, in eq. (B5) 
^2 


FCN difference function for iterative solution for TOUT 
TPRE previous value of TOUT 


Program listir^ . - The FORTRAN listings for function subprograms CP, H, T2H, 
S, T2S, and HF and for subroutine THERMO are as follows; 

PUNCT ICN CP(T,F ) 

CALL TLEUMO 
CP = X 
RETURN 
END 
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o o n o r> o o 


FUNCTION H(TltT2fF) 

CALL thermo (Tl,T2fF,X,2) 
H=X 

RETURN 

ENO 


FUNCTION T2H(TliF,X) 

CALL THERMO ( T I ,T 2 f F , X , 3 > 

T2H=T2 

RETURN 

END 


FUNCTION S(Tl,T2fF) 

CALL THERMO ( T 1 ,T 2 » F , X , 
S=X 

RETURN 

END 


FUNCTION T2S(Tl,F,X) 

CALL THERMO ( T I ,T 2, F, X , 5) 

T2S=T2 

RETURN 

END 


FUNCTION HF(T1,T2) 

CALL THERMO ( T 1 , T 2 , 1. f X , 6 ) 
HF=X' 

RFTimN 

END 


SUBROUTINE THERMO ( TI N , TOUT F OA , 

THERMODYNAMIC PROPERTY SUBROUTINE 

IND r.IVEN 

1 
2 


FUNC, IND) 


CALCULATED 


T IN,FOA 
TIN, TOUT, FOA 


C.P(TIN,FOA) 

H(TOUT,TIN,FOA) 



3 


C 

C 

c 

c 

c 

c 

c 


TINfFPa ,H(T(KJT, TIN,FQA) 
1 IM,TOUT,Fnfl 


5 TIN,FnA,S(TmiT,TIN,FnA) 

6 lIN,Tn'JT 
COMMnN FOC, TTOL, W 


TOUT 

S(TC JT,TIN,FOA) 
TOUT 

HF ( TOUT , T I N 1 


THEPMOCYNAMi: PROPERTY FiJNCTIONS 


CP( T, F )= ( l.0/( 1.0+F+W) ) 

1 *( A + B*T+C*T’f'*2<-0»T«*3«- E*T>!'*4<-( F/ ( 1 .O+HOC ) I 

l#( Al* ( 4r2+B02*T+C02*T’»*2«-002*T*=!!3«-En2*T«^4 ) *■ A2 ♦ ( AC 02 +BC02 *T+CC02* 
2T**2+DCL'2*'T>!'* 3<-FCD2*T**4J 
1) ♦( A3*F/( l.O + HOC »+W) 

1 AH20<-RH2n*T4-CH2n<‘T**2+[)H20*T*#3 + 

3EH20*T**4 ) ) 

H(T2,T1,F)= ( 1.0/( 1.0 + F + wn 

1 * (A*( T2-T1) +8*(T2«’<'2-Tl<'*2)/2.0 *C*(T2 

1#*3-T l**3)/3.0 +P*( T2**4-T1’)'’«'4) /4.0 4-6* ( T2**5-T l=f*5 ) / 5 . 0 +(F/(l.O + 
2H0C ) »* ( A 1*(A0 2* { T2-T1 > 4-B02# ( T2**2-T 1* *2 ) /2 . 04 C 02« ( T 2**3-Tl **3 ) / 3.0 
34nn2»( T2*>i‘4-T 1#*4 )/4.04t02«( T2*’«''5-T1«=«5 ) /5.0) 4 A 2^ ( AC02* (T2-T 1 ) 4 
4PC0 2* ( T2«*2- Tl*#2)/2.04CC02*(T2-i'>i«3-Tl*’*'3) /3.0 4DCO2’f'(T2#*4-Tl**4)/ 
54.04EC02*( T2t'*5-Tl*X'S ) /S. 0) 

U 4( A3*F/{ 1 .04HQC ) 4W) 

5 *(AH20><=(T2-T1)4BH20*(T2**2-T 1«*2)/2. 

64C H 20 *( T2** 3- Tl** 3) /3 . 040H20* ( T2 4-T ) /4. 04 EH20 *( T 2**5-T 1**5 ) / 

75.0 n 

S(T2,T1,F)= ( 1.0/( 1.04F4W) ) 

1 *(A*ALOG( T2/T1 )4B*{ T2-T1 ) 4C*( T2**2-T 1**2 

1)/2.04C*( T2**3-T1**3» /3.04E*( T2**4-T1**4) /4.0 4 (F/ ( 1.04H0C ) )*( Al* 
2( AD2*AL0G(T2/T1|4B0 2* ( T2-T1 ) 4C02* ( T2**2-T 1**2 )/ 2.0 4 002 *(T2 **3-T 1** 
33 >/3.04L02*( T2**4- Tl**4) /4. 0) 4 A2* ( AC 02 * ALOG ( T2/T1) 4BCG2*( T2--T 1 > 
44CC02*(T2**2-T1**2) /2 . 04DC02*( T2**3-T 1**3 ) /3. 04EC02 *( T 2 **4-T 1**4 
5 )/4.0) 

1 ) 4( A3*F/( 1 .04H0C I 4W) 

5 *(AH20*AL0G( T2/T1)4-8H20*(T2-T1)4CH20*(T2**2-T 1**2) / 

62.04DH20*(T2**3-T1**3) / 3. 04EH20* ( T2**4-T 1 **4 ) /4.0) ) 

HF(T2,Tl) = ( 1.0/( 1.04H0Cn* 

1 (A1*«A02*(T2-T1 )4B02* (T2**2-T1«*2) /2. 04 C02* < T 2**3-Tl **3 ) / 3.0 

24[)G2*( T2**4-T 1**4 I/4.04F0 2* f T2**5-T1**5) /5.0) 4 A2* ( ACG2* (T2-T1) 4 
3BC0 2*( T2**2-T1**2 )/2.04CC02*(T2**3-T1**3) /3 .0 40C02 * (T 2 **4-T 1**4 ) / 

44.0 4EC0 2*(T2**5-T1**5) / 5. O ) 4 A3* ( AH2 0* ( T2 -T1 ) 4 BH 20* ( T2 **2-T 1**2 )/2. 
54CH20* (T2**3-T 1**3) /3 . 040 H20* ( T2**4-Tl**4 ) /4. O 4 EH2 0*( T 2**5-T 1**5 ) / 


THERM - EFN SOURCE STATEMENT - IFN(S) 


65. 0 ) ) 

A = .2AC62 
B = -.C 17724E-3 
C = .C28C56E-6 
0 = - .0 12662F-9 
E = .CC12012E-12 
A02 = .20334 
B02 =.C2968E-3 
C02 =.CC3997lE-6 
002 = -.GC53B42E-9 
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PVl = .C0G7676''+E-12 

4C02 =.11C97 

«.C72 =.2lllOE-2 

CCQ2 =-.0881^0E-6 

Dr.02 =.C18003E-9 

ECQ2 =- .00 143 I7r- 12 

AH20 =.44266 

RH20 =-.033168?-3 

CH20 =.C8776lf-6 

0H2n =-.C24552E-9 

EH20 =.C021734E-12 

A1 =- 2.66 44*( l.C+2. 9787*HOC I 

42 = 3.6644 

A3 = 8.9 365#HnC 

Gn T3 (10, 20, 30,40, 50,60 , IND 
10 FIJ\JC = CP( TIN,F0A ) 

RETUPN 

20 FlJNC=H(TnUT,TIN,FGA) 

RETUR N 

30 Tni)T = T IN-*-FIJNC/CP(TIN,F0A( 

31 FCN=FUNC-H( TOUT, TIN ,FGA) 

DFCN=-CP( TaHT,FOA ) 

TPRE=TnUT 

TOUT=TCLT-FCN /OFCN 

I F{ A8S(TnuT>TPRE) .GT.TTOU GO TO 31 
RETURN 

40 FIJNC=S( TOUT, TIN ,FOA ) 

RETURN 

50 TOUT=T IN*EXP( FUNC /C P ( TI N, FOA ) ) 

51 FCN=FUNC-S(TnUT,TIN,FOA» 
OFCN=-CP(TOUT,FOA )/TGUT 
TPRE=TOUT 

TnUT=TruT-FCN70FCM 

IF(ABS(T(3UT“TPP.e) .GT.TTOL) GO TO 51 
RETURN 

60 FUNC=HF( tout, TIN) 

RETURN 

END 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, October 3, 1974, 

501-24. 



APPENDIX A 


Ki,K2,K3 


r 

SFC 

T 

w 


SYMBOI5 

constant in heat capacity equation, eq. (B4), J/(kg)(K) (Btu/(lb)(°R)) 
constant in heat capacity equation, eq. (B4), j/(kg)(K^) (Btu/(lb)(°R^)) 

constant in heat capacity equation, eq. (B4), J/(kg)(K^) (Btu/(lb)(°R^)) 

heat capacity, J/(kg)(K) (Btu/(lb)(°R)) 

constant in heat capacity equation, eq. (B4), J/(1^)(K^) (Btu/(lb)(°R^)) 

constant in heat capacity equation, eq. (B4), j/(kg)(K^) (Btu/(lb)(°R^)) 

fuel-air ratio 

enthalpy function Ah as defined by eq. (B7), J/kg (Btu/lb) 
enthalpy function Ah^ as defined by eq. (C4), J/iig (Btu/Ib) 
enthalpy difference or specific work, J/kg (Btu/lb) 
heat of combustion of fuel, J/1^ (Btu/lb) 

enthalpy correction function as defined by eq. (BIO), J/kg (Btu/lb) 

dimensional constant, 1 (778 ft-lb/Btu) 

ratio of turbine power to compressor power 

constants in heat capacity equation, eq. (B5) 

molecular weight 

ambient air absolute humidity 

number of fuel compressors with intervening intercoolers 

number of shafts 

power, kW (Btu/sec) 

absolute pressure, N/cm (psia) 

gas constant, J/(kg)(K) (ft)(lbf)/(lbm)(°R)) 

ratio of component exit to inlet total pressure 

specific fuel consumption, kg/(hr)(kW) (lb/(hr)(hp)) 

absolute temperature, K (°R) 

mass flow rate per pound of inlet air, 1^/sec (Ib/sec) 


46 



X atom ratio of hydrogen to carbon in fuel 

y mass ratio of hydrogen to carbon in fuel 

7j efficiency or effectiveness 

$ entropy function ^<p as defined by eq. (B12), J/(kg)(K) (Btu/(lb)(°R)) 

entropy function t^cp^ as defined by eq. (C8), J/(kg)(K) (Btu/(lb)(°R)) 
A(^ constant-pressure entropy change, j/(kg){K) (Btu.'(io)(°R)) 

Subscripts: 
air air 

B burner 

C compressor 

CO 2 carbon dioxide 

c turbine coolant 

cold cold side 

cv conversion 

cy cycle 

D diffuser 

est estimated 

ex exit 

f fuel 

HgO water vapor 

hot hot side 

I inlet component 

INT intercooler 

id ideal 

in inlet 

j shaft j 

1 recuperator leakage 

m mixed 

net net 

0 2 oxygen 



o overall 

p polytropic 

R recuperator 

r reference 

sh shaft 

T turbine 

tot total 

0 ambient condition, see fig. 5 

1 initial; station between inlet and compressor, see fig. 5 

2 final; station between compressor and recuperator, see fig. 5 

3 station between recuperator and primary burner, see fig. 5 

4 station between primary burner and turbine, see fig. 5 

5 station between turbine and diffuser, see fig. 5 

6 station between diffuser and recuperator, see fig. 5 

7 station at recuperator hot-side through-flow exit, see fig. 5 

8 station at recuperator hot-side exit after leakage mixing, see fig. 5 

Superscipt: 

' absolute total state 
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APPENDIX B 


DERIVATION OF COMBUSTION -GAS PROPERTY EQUATIONS 


Reaction Stoichiometry 


The fuel used has the composition CH^, and the general combustion equation with 
associated formula weights is 


CHx 


r~ ^ 

12.010 + 1.008 X 




(Bl) 


The variable used in the computer program to ejq)ress fuel composition is the mass 
ratio y of hydrogen to carbon. In terms of this mass ratio, the atom ratio x is 


^ 12.010 y ^ 11.914683 y (B2) 

1.008 


and the formula weights associated with the components of equation (Bl) are 

CHjj = 12.010(1 + y), ©2 = 32.000 + 95.317460 y, COg = 44.010, and HgO = 107. 327460 y. 

If one mass unit of air, with an associated m mass units of water vapor, is reacted 
with f mass units of fuel, the reactant and product masses w per mass unit of air are 
as shown in table IV. The composition of the combustion gas used in the cycle analysis 
is defined by the product masses shown in table IV. 


Heat Capacity 

The heat capacity of the combustion gas is 






p ^components 
1 + f + m 


(B3) 


The heat capacity of each component is expressed in the form 

c = A + BT + CT^ + DT^ + ET*^ (B4) 

The coefficients A, B, C, D, and E were obtained by a least-squares regression analy- 
sis of the tabulated data of reference 1. These coefficients and the molecular weight of 
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each component are given in table V. 

Equation (B3) for combustion gas heat capacity is evaluated by substituting for w 
according to table IV and for Cp accordii^ to equation .(B4). Performing. the siunmation 
yields 




+ y/ 




+ Ko A 


2\"CO. 


+ B- 


CO, 


T + C 


CO, 


+ 


D 


CO, 



+ 







(B5) 


where 

A^ir = 1.00821 J/(kg)(K) (0. 24062 Btu/(lb)(‘^R)) 

Bair = -1.33 6 75X10"^ j/(kg)(K^) (-0.017724x10"^ Btu/(lb)(°R2)) 
C^. = 5, 16637x10"^ J/(kg)(K^) (0.038056x10"® Btu/(lb)(°R®)) 

Aiir 

= -3.09412x10"'^ J/(kg)(K^) (-0.012662x10"® Btu/(lb)(°R‘^)) 

drlP 

= 5.72336x10"^^ J/(kg)(K®) (0.0013012x10"^^ Btu/(lb)(°R®)) 

aIP 

Kj = -2.6644 (1 + 2.9787 y) 


Kg = 3.6644 


, Kg = 8. 9365 y 

The Og, COg, and HgO coefficients are as listed in table V. The coefficients sub- 
scripted air were obtained by combining the component coefficients accordii^ to the 
air composition shown in table IV. 
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Enthalpy Change 


The change in gas enthalpy resulting from a change in gas temperature from Tj to 
T2 is 


Ah = 



(B6) 


Substituting equation (B5) into equation (B6) and integrating yield 



For brevity in the cycle analysis procedure, equation (B7) will be expressed as 

Ah = H(Tj,T2,f) (B8) 

In order to facilitate burner calculations, it is desired to extract the fuel-air ratio 
f from the function H. To do this, equation (B8) is written 
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Ah = 1 (1 + ni)H(T.,Tp,0) + fI(Tj,T 2 ) 

1 + f + m ^ ^ ^ ^ 


Comparison of equation (B9) with equation (B7) shows that the function I(Tj, Tg) is 


KTpTg)^ 


1 + y. 


Ao 2^'^2 - ^ 


Bn 

^2 L2 r^2\ . ^2 


(t|-t2) + -^(t|-t5) 


" - Tt) - (i - i) + Kj Apo^(T2 - Tj) . (t2 - 


3 V ^ V 4 


-(T|-Ty + -^(T|-T5) 


♦ ^3 Ah^q(T2 - Tj) > (t| - T?). (t| - t5) 


(t«-t«)*_^(t|-t5) 


(BIO) 


Constant-Pressure Entropy Change 

The change in gas entropy resulting from a change in gas temperature from to 
T 2 at constant pressure is 



(Bll) 


This function is used in the evaluation of ideal (isentropic) and polytropic processes in 
turbomachines. Substituting equation (B5) into equation (Bll) and integrating yield 
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For brevity in the cycle analysis procedure, equation 

A<p = #(TpT2,f) 


Tj) + 1 - 

^ 2 




(B12) will be expressed as 


(B12) 


(B13) 


Molecular Weight 


The molecular weight of the gas is equal to the weight of the gas divided by the total 
number of moles (sum of the moles of the components). Therefore, molecular weight 
can be expressed as 


M = 


1 + f + m 

y(~) 

components 


(B14) 
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With w obtained from table IV and M from table V, equation (B14) expands to 


1 + f + m 

0.034522 + 0. 24802 y^— ) + --- - - 

ll + y/ 18.016 
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APPENDIX C 


FUEL PROPERTIES AND COMPRESSION 
Thermodynamic Properties 


The heat capacity of the CH^ fuel is expressed as 


^p,f - ®f"^ + 


(Cl) 


The coefficients A^, B^, and Cj depend on the particular fuel being used and must be 
provided as program input. 

The change in fuel enthalpy with temperature change, in accordance with equation 
(B6), is 



(C2) 


In accordance with equation (Bll), the change in fuel entropy resulting from a tempera- 
ture change at constant pressure is 


'^2 

= Aj In — + BjCTg 
”^1 




(C3) 


In a similar manner as was done for the combustion gas in appendix B, equations (C2) 
and (C3) will be expressed 


Ahj ^ Hf(Tp Tg) (C4) 

and 

A<p^ = $j(Tj,T2) (C5) 


Fuel Compression 

In some cases it may be necessary to compress the fuel in order to get it into the 
burner. For such a case, it may be desirable to chaise the cycle output for the fuel- 
compression power. The program provides for the computation of fuel compression 
work with or without intercooling. Although it is recognized that a gaseous fuel may 
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deviate significantly from ideal-gas behavior, an ideal-gas analysis is used because the 
fuel-compression power is small as compared with the cycle output power. 

The fuel inlet temperature fuel inlet pressure pj fuel compressor poly- 

tropic efficiency i]q j p, and number of compressors n^ ^ with intervening intercool- 
ers are specified by the program input. It is assumed that the fuel must be compressed 
to a pressure defined by the cycle pressure ratio r 


cy 


Pf,C,ex = P0^cy 


(C6) 


Each fuel compressor is assiuned to have the same pressure ratio 


/Pf.C.exV /^C.f 

\ Pf,in / 

and each intercooler is assumed to cool the fuel back to the fuel inlet temperature 
Tj The fuel-compressor exit temperature Tj ^ is then obtained from 


(C7) 


♦fKln’T;,c,ex) 


1 ^f 

^ i In rp f 

J 


(C8) 


'^C,f,p 


and the fuel-compression specific work is 

C " "C, f”f('^f, in’ "^f, C, ex) 


(C9) 
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TABLE n. - OUTPUT FOR ILLUSTRATIVE EXAMPLE 


uCvtr SV5 .:l 1 ^ r, AT A.'-UcNT C UNf. r TI ONS ZSfl.OO nfG'KcLVIN AMO 10.130 N/CM2 

iNLt.r TOTAL IFi-SSO-l F'iC" 1 V_: ^ v = j . 900 » T‘J.?dI'Ic. tX^A'»ST HICPUSE? TOTAL PRESSURt « cCOV £P Y = 0. 990i EXIT STAT/TOT PRESS RATIf)=0.980 
HtCUP£y ATOl; ;;fFc( SS = .).O00, O^L'.) S I L‘ c TOTM. u-^cSS»ip.t o t-CnVtT Y = 0.930* HTT SIDE TOTAL PRESS'JRf- » cCOVfcR Y= 0. 96 0i BYPASS PLOWED. 

PULL Hi AT inn VAl.U- = ‘Vij‘'JJCJ. AT 208. D-', K-zLVl". , H/C 9 AT 1 1>=0. 30770 
ruit. OaS f.P= ft69.'.u00 *■ ♦•00 <■ 7 l.352Si-03 * T**2 

FUlI. M(L wnT= Ift./i-.o f.,i:-L IUL:.V Tf.-'-.ps 219. » FUt:L I'JLrT C SS I CN= -35.0# TUEL COMPRESSION fepP. *0.880 

SHAFT /! -HAi A ft ? S Si U- S THAT SUPPLY u.6jJu UP TOTAL COMPRESSION AND 2 TlJPRINtS 

Covucusjf^ R PPOVIO.^S 0.3000 "tP SHAFT CCMPR*SSION WITH fePF I CIcNCVaO . 880 

FLOW rUTLf- If-'-. C('M'‘k.‘iSSiU 3 IS CUULttD TO J06. Opr, KrLViN WITH TOTAL PRESSURE R cCOV EP Y =0 .9 80 
CruFPi.SSCF 3 PFOVlOt'S 0.2*500 HE SHA= T C»lMPkLSSlUN WITH dFFI CICNCY=0.flfl0 

FLOW rNTLFIM'*. COMPFESS'JW 2 IS C.5JL:;U TO 306. ULG KELVIN WITH TOTAL PRESSURE R ECOVFR Y* C .980 
r'jMPRFSSCJC 2 PfOVIUcS O.2S00 liF SHAFT COMoPcSSln-.l WITH EFFI CIdNCY = 0. 380 

FL J>i LNTcFINf. r.(]MpK.-S>OP 1 IS CiiLLcO TU 306. DEC KELVIN WITH TOTAL PRESSURE RFCnVERY = C.980 
Cl'IPKtSSnP I PfnVIDcS C.2000 of shaft COMPEtSSION WlTH tFF|CIENCV*0.880 

FLOW FNTcr- INS TUPhlUr I IS HKATtO TO 1 389.0 DEO KELVIN WITH' COMBUST ION EFFICIENCY =C.980 AND TOTAL PRESSURE RECOVER Y*0. 970 
TURHiUc I FpowindS O.AOiiO '^F SHAFT cYdANSION wITH E FFI C I ENC Y=u. 900 
TiJf^HlK'c 1 IS COnLED WITH J.u25o OF CnMPktSSOR EXIT FLOW 

FLOW cMEFIMr. TliKHlNE 2 IS HEATlO T(l 1389.0 DEG K.-LVIN WITH C0«8lSriDN EFFICIENCY *C.9pO AND TOTAL PRESSURE R£C0VERY«0.970 
TURIU.'-E 2 PKOVIUhS O.oJJJ OF SHAFT EXPANSION WITH fPFK ! ENCY=0.900 
TijF'iiNi 2 IS rur.Lrn with j.0250 cf cqmprlssor exit flow 

SHAFT 1 HAS 1 COYPR r.SSnPS THAT SUPPLY C.AOOU OF TOTAL TOMPRcSSION AND 2 TURBINrS THAT DELIVER. 1.0000 • COMPRESSOR SHAFT POWER 

FLOW cUTEFiMr, (.Oms^LSSOR 1 IS CJULEO TO 306. 'OcG KELVIN WITH TOTAL PFESSURE P. ECOVtR Y * 0.9 80 
COMPCr-SSr--. I PrnvinHS UOOOO of shaft CCMorlsSION with cFfi cI£NCy»0.880 

FLOW tNlEhlUG TUPPlUe 1 IS HfATED TU l389.0^3cG KELVIN WITH COMBUSTION EFFICIENCY »C.980 ANO TOTAL PRESSURE ReC0VERY*0. 970 
TURMlUb 1 PfOVIOcS 0,r*000 OF SHAFT PCWEH W I T H t FFI C I cNC Y*0.90U 
TURBINE 1 !S f fiOLl-0 WI Ttl J.0250 CF COMPRESSOR EXIT FLOW 

Finu i:N7rRINfi TUBpIUF 2 I .S hlaTcD Ti! 1389.0 OL-G KtLVIN WITH COMailSTION EFFICIENCV »C.980 AND TOTAL PRESSURE R6COVERYs0.970 
TUFBlNr < OfOVIOcS O.SOJO OF SHAFT POWcR WITH CFFI C I ENC Y*0.900 
TURF INI 2 IS rnOLtO WITH 0.029. CF COMPRcSSOP. EXIT FLOW 


OUTPUT P.JWi-F 

EnilALS 

0.9500 HF 

SHAFT POWER 


TUR RCMACHINFPY 

EFFlClENCItS ARE PDLYTRDPI 

PPLSS. 

N.tT 

Ni:T 

FUJI 

CYC LE 

F!JEL 

POW::R VALUES REFER 

RAl iO 

pr.wt ^ 

PL* '4r:9 

c wsuMo. 

EFF 

POWt.R 

TU ONE KO PcR SEC 


K'v 

HP 

K0/f3«/Kw 


KW 

OF inlet AIR FLOW 

5.0J 

255. P6 

293. 3J 

J.2209 

0.377O 

0. 


10. 0) 

2FF.26 

520. 95 

0.1898 

0.9379 

0. 


15 ,U0 

959. 79 

616. 87 

0.L819 

0.958?. 

0. 


20.0U 

507.87 

691. 99 

0.1779 

0.9689 

0. 


25.00 

592.73 

729. 56 

U.1751 

0.9 796 

0. 


30.00 

572.10 

767. 62 

0. 1 736 

0,9786 

0. 


35.00 

595.9,3 

798. S3 

0. 1 726 

0.9819 • 

0. 


90.00 

a -5. 19 

825. 91 

0.1 719 

0.9839 

0.0031 


95 .u-J 

<22.20 

398. 27 

0.1719 

0.9899 

0 .0060 


50.00 

<9h.97 

868. C8 

0. 1710 

0.9860 

0.1919 


55.00 

<59.85 

88 5. 3 7 

0. 1 707 

0.9867 

0.6053 


60. Ou 

<71.90 

9U0. 87 

0. 170.5 

0.9873 

0.9981 


65.00 

Cf 1.85 

919. 89 

C. 1 7 09 

0.9877 

1 .381 1 


70.00 

<91.37 

927, 69 

a , 1 703 

0.93BO 

1.7971 


75. JO 

7CC.10 

939. 33 

C. 1702 

0.9882 

2.1000 


PO .00 

7C5.19 

950. 19 

V.1702 

0.9883 

2.9930 


85,00 

? 1 5 . p 2 

96C. 20 

C. 1 701 

0.9689 

2-7726 


90.00 

722 .59 

999. 52 

C. 1701 

0.9889 

3.0955 


95.00 

7 2 9, Jo 

973. 23 

C. 1701 

0.9889 

3.9059 -• 


UO.Oj 

7?5 . 19 

989. 39 

0.1702 

0.9889 

3.7105 




TABLE IV. - COMBUSTION REACTION STOICHIOMETRY 


Component 

Reactant 
mass per 
unit mass 
of air 

Mass formed (+) or consumed (-) 
by reaction per unit mass of air 

Product mass per unit mass of air 

^2 

0.7552 

0 

0. 7552 

n. 

.2314 

32.000 + 95.317460 y , 

n OQi/i 2. 664446 + 7. 936508 y f 






. 12.010(1 + y) 

1 + f 

A 

.0129 

0 

0.0129 

CO, 

.0005 

1 44.010 j 

0.0005 + ^ ’ ^ f 



12.010(1 +y) 

1 + y 

H,0 

m 

107.327460 y , 

8.936508 y , 
m + — I 



12.010(1 + y) 

1 + y 

CHx 

f 

-f 

0 

Total 

1 + m + f 

0 

1 + m + f 


TABLE V. - COEFFICIENTS AND MOLECULAR WEIGHTS OF COMPONENTS 


Compo- 

nents 

Coefficients in equation (B4), SI units (U. S. customary units) 

Molecular 

weight, 

M 

A 

B 

C 

D 

E 

N, 

1064.69 

-246. 550x10'^ 

647.357x10'® 

-365. 884x10'® 

65.4721x10'^^ 

28.016 


(0. 25410) 

(-0.032690x10'^) 

(0.047685x10'®) 

(-0.014973x10'®) 

(0.0014885x10'^^) 


°2 

852.001 

223. 848x10'^ 

122.142x10'® 

-143. 788x10'® 

33, 7649X10'^^ 

32. 000 


(0. 20334) 

(0.029680x10'^) 

(0.0089971x10'®) 

(-0.0058842x10'®) 

(0.00076764x10'^^) 


CO2 

464.968 

1592. 13x10'^ 

-1196.56x10'® 

439.926x10'® 

-62.9737x10'^^ 

44.010 


(0. 11097) 

(0.21110x10'^) 

(-0.088140x10'®) 

(0.018003x10'®) 

(-0.0014317x10'^®) 


H,0 

1854. 76 

-250.057x10'^ 

1191.42x10'® 

-599.959x10'® 

95.5976x10'^® 

18. 016 


(0.44266) 

(-0.033155x10'^) 

(0.087761x10'®) 

(-0.024552x10'®) 

(0. 0021734x10'^®) 


A 

521.156 

0 

0 

0 

0 

39. 944 


(0. 12438) 

(0) 

(0) 

(0) 

(0) 



61 




























High-pressure 

compressor—. 


Combustor - 


- High-pressure 
I turbine 


^ / 
Fuel-' 

' 

Intercooler 

Reheater 




Hot gas 


exhaust 


Low-pressure ' 
turbine 


Figure 1. - Open-cycle system with Intercooling and reheating. 
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Figure 2 . - Thermodynamic diagram for open-cycle system with 
intercooling and reheating. 
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Figure 3. - Effect of intercooling and reheating on cycle efficiency 
at optimum pressure ratio. Compressor polytropic efficiency, 

88 percent; turbine polytropic efficiency, 90 percent. 
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(a) Maximum number of components. 



(b) Minimum number of components. 
Figure 5. - Open-cycle power system schematic. 
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